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ABSTRACT 



The power spectral tcchnlquca of generalized harmonic analysis 
were applied to determine the phase and amplitude of the transfer 
function between the vertical component of turbulence and the 
resultant pitch rate of a small swept wing jet fighter. The power 
spectrum analysis was based on actual flight test data. In addition 
the amplitude and phase of the aircraft's transfer function are 
determined theoretically by^..two different methods and compared with 
the results of the power spectrum analysis. It was established that 
the power spectrum method, as applied to random turbulence, can be 
used to determine an aircraft's response to this random input on a 
continuous basis and that power spectrum techniques have considerable 
potential value in the study of aircraft gust loading. 
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INTRODUCTION 



Atmospheric turbulence and the resulting aircraft disturbances 
from stabilized flight due to this turbulence are random and thus 
aperiodic and non-dlsslpatlve. These disturbances can only be de> 
fined on a non-preclse basis with certain probability distributions 
and average values since they are random in nature. This leads to 
the application of probability and statistical concepts, which. In 
combination with classical Fourier techniques, form a general 
approach for the analysis of problems of this type. This technique 
Is termed the power spectrum method. 

Power spectral methods were derived from economic time studies 
at the beginning of this century. The initial applications were also 
directed toward turbulence studies and the Brownian motion of parti- 
cles. During World War II power spectral methods were used to analyze 
problems In servomechanisms and communications; however, the applica- 
tion of power spectral analysis to the field of aeronautics, and in 
particular to the study of the dynamic response of aircraft to tur- 
bulence, Is relatively new and still In the early stages of development. 

Other recent studies In the field of aeronautics Involving 
power spectral techniques have been applied to: 

a. Wind tunnel turbulence 

b. Gust loads on aircraft 

c. Acoustic Noise 

d. Buffeting loads 

e. Electronic noise in guidance aAd gunnery systems 

In particular, considering the gust loading and dynamic response 
of aircraft, statistical Inference Is a realistic approach since the 
problem is essentially one of a random process phenomena. Perhaps a 
more meaningful presentation of the effectiveness of the power spec- 
trum method can be made If the earlier discrete gust design methods are 
mentioned in comparison. 
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The earliest discrete gust approach involved the assumptions of 
a rigid airplane f a ramp type gust, and a specified maximum vertical 
velocity. Subsequent discrete gust design methods were adjusted for 
different wing chords and pitch characteristics between different air- 
craft by expressing data in terms of the mass ratio . The dis- 
crete gust concept was based primarily on the premise that the rela- 
tive loads for single isolated gusts are a measure of the relative 
loads in a sequence of gusts. 

However, in recent years a more general approach to turbulence 
analysis has been needed in order to account for the greater vstrieijt*' 
and degree of aircraft flexibility, mass distribution, configuration 
and dynamic response characteristics, as well as the varied operating 
speeds and altitudes. Therefore, the power spectral techniques of 
generalized harmonic analysis were developed for application to tur- 
bulence and dynamic response studies. 

As conditions changed and the variety of aircraft and operating 
altitudes/airspeeds increased, the concept of the discrete gust design 
method was broadened to include flexible bodies as well as rigid 
bodies. Although the flexible aircraft was mathematically flown 
through discrete gusts of varying wavelengths, the discrete gust con- 
cept did not adequately account for the continuous nature of atmosf- 
pheric turbulence, particularly as applied to lightly dasiped dynamic 
systems. Latest design philosophy embodies both the discrete gust 
concept and the power spectral approach, with the design initially 
analyzed and checked using the discrete gust method and then examined 
in greater detail using the power spectrum analysis. Further develop- 
ment of the power spectral method may result in its use in the future 
entirely for atmospheric turbulence studies. 

The main advantages of the power spectral approach to turbu- 
lence analysis are outlined in Ref. 1: 

1. It provides a sure realistic representation of the “ ' 
continuous nature of atmospheric turbulence. 

2. It allows airplane configurations and response character- 
istics to be taken into account in a rational manner. 

3. It allows more rational consideration of design and 
operational variations such as configuration changes, mission 
changes and airplane degrees of freedom. 
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Most of the present day applications of power spectral methods 
are concerned with the analysis of the magnitude and distribution of 
the random data; for example, the intensity of atmospheric turbulence 
and its variation with frequency. However this thesis considers a 
more unique application of power spectral techniques by establishing 
the transfer function between pitch rate and turbulence for the 
F-IOOC airplane based upon the analysis of a random input, w , the 
vertical component of atmospheric turbulence. This application goes 
one step beyond the analysis of random data by considering its effect 
upon the dynamic response of an airframe. The particular transfer 

A 

function, ~ , is considered to provide a basic indication of the 
F-lOO airplane's dynamic response characteristics in a turbulent en- 
vironment. To measure the accuracy of the results obtained through 
the power spectral analysis, a comparison between , as deter- 

mined from this method, is made with the values of jL. , as deter- 

‘^9 

mined by solving the longitudinal equations of dynamic motion for 
the F-IOOC airplane. 
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DISCUSSION OF THE CHARACTERISTICS AND ANALYSIS OF RANDOM DATA 



Any recorded data which describe a physical process can be 
broadly classified into either one of two categories: deterministic or 

non-deterministic. Deterministic data are information which describe a 
process which can be represented and predicted with reasonable accuracy 
by exact mathematical relationships. Some common examples of this type 
of process are the electrical potential across a condenser as it dis- 
charges through a resistor or the motion of a satellite in a known 
orbit around the earth. Examples of non-deterministic data are the 
variation of an aircraft's vertical velocity in turbulence and the 
variation in the magnitude of the lift developed on an airfoil encounter- 
ing buffeting in rough air. Due to the inherent randomness of the 
turbulence this non-deterministic data cannot be described, nor can its 
behavior be predicted by explicit mathematical relationships. In 
practical terms, the decision as to whether or not physical data are 
classified as being deterministic or non-deterministic (random) is 
based upon the ability to reproduce the data repeatedly through con- 
trolled experiments. If the experiment or phenomena which produced 
the data can be repeated with the same results, within the tolerances 
of experimental error, then the data obtained are classified as deter- 
ministic. Conversely, if identical results cannot be obtained upon 
repeating an experiment, then the data are considered to be non- 
deterministic or random in nature. 

In considering random data, therefore, it is expected that an 
particular observation producing this data will effectively represent 
only one of many possible outcomes which might have occured as a result 
of this experiment. For example, assume that the lift developed on an 
airfoil during buffeting in flight is recorded as a function of time. A 
specific lift time history will be obtained as shown in Fig. 1. 

However, if a second airfoil of identical shape is flown through a 
similar type of disturbance simultaneously, it is likely that a different 
lift time history would result. Therefore, the lift time history for 
any one airfoil encountering a particular finite disturbance is only 
one example of an infinitely large number of time histories which might 
have occurred. If the random buffeting is observed over a finite time 
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L/f^r 













Figure 1. Sample records showing the lift developed on an airfoil 
in turbulence. 
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interval and is graphically or mathematically represented by a single 
time history, then this time history is referred to as a sample record 
by Bendat and Piersol (2). The collection of all possible sample 
records of lift versus time which the random buffeting might have pro- 
duced is called a random process or stochastic process in Ref. 2. 

As a result, a sample record of observed data for a random physical 
occurrence can be considered to be only one physical outcome of a 
random process. 

The properties of a physical experiment or occurrence which 
are described by a random process can be determined from the average 
values of the experimental data as computed over the range of sample 
records obtained. A collection of sample functions which comprise a 
random process is illustrated in Fig. 2. By taking the instantaneous 
values of each of the sample functions at time t^, summing the values 
and dividing by the number of sample functions, the mean value of the 
random process for a particular time t^ may be computed. In addition, 
a Joint moment or correlation between the values of the random pro- 
cess at two different times can be computed by taking the average of 
the product of the instantaneous values at two different times t^ and 
t^ + Z/ . The product obtained in the manner described above is termed 
the auto-correlation function. The mean value, auto- 
correlation function/^ defined as: 

V 

^ /V->co ^ 



•"<* /?^ (t,^ t, ^7) ii Z X, (t) -X. (t, ^7) m 

A/-*oo ^ 

According to Ref. 2, the random process Is stationary if 

do not vary as time t^^ is varied. In a 
more general sense, when the properties of a random process, as com- 
puted over short intervals, do not vary significantly from one interval 
to the next, then the process is considered to be stationary. 
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Figure 2. 



Applying the concept of stationarity to atmospheric turbulence. 

It has been shown In Ref. 1 that the conditions of stationarity and 
homogeneity apply to large scale patterns of air motion within regions 
of one hundred miles and time durations of one hour. In cases of severe 
turbulence, such as encountered in a thunderstorm, these time and dis- 
tance estimates would not be valid. However, the tlsie duration for 
homogeneous and stationary conditions may still be almost five minutes 
in length. It Is usually possible to determine the properties of a 
stationary random process by computing time averages over specified 
sample functions In the ensembles (collection of sample functions). 

Consider the sample for the random process Illustrated In 
Fig. 2. The mean value, , and the autocorrelation function, 

y ( -^ ) » for sample function, according to Ref. 2, are: 






(3) 






/^x.(t)-xj(t^y)dt ( 4 ) 

Random processes may be classified as shown in Fig, 3, repro-* 
duced from Ref, 2, as being either stationary or non-s6at ionary , and 
stationary random processes may be further subdivided as either 
ergodic or non-ergodic. The random process is ergodlc if it is sta- 
tionary and an d/^j^(Z4)do not vary when computed over 

different sample functions. Therefore, all properties of ergodic ran- 
dom processes can be determined by computing over a single sample 
function. In actual physical phenomena, as outlined in Ref, 2, most 
random data, including atmospheric turbulence, are ergodic so that 
the properties of the data can be computed from time averages of 
individual sample records. 

Two types of statistical functions are used to describe the 
basic properties of random data; (a) the autocorrelation function 
and (b) the power spectral density function. The autocorrelation 
function and the power spectral density function provide an analysis 
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Figure 3. Classifications of random data (Reproduced from Ref. 2.) 
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of the amplitude of the random data with respect to tine and frequency 
respectively. These functions are Fourier transform pairs, and there* 
fore provide an analysis of the data in two different forstats. 

The autocorrelation function for random data provides an indi* 
cation of the general dependence of the amplitude of the data at a 
particular time on the value of the data at a different time. An 
autocorrelation measurement can be used for detecting deterministic 
data i^ich nuay be hidden or obscured in a random background. 

As an extension of the autocorrelation function, the power 
spectral density function describes the frequency composition of the 
data in terms of the spectral density of its mean square value. 

In Chapter IV an explicit and more detailed presentation is 
included which outlines both the analytical and digital computer steps 
involved in the power spectrum analysis as used to determine the pitch 
response of the F-IOOC airplane to atmospheric turbulence. This 
analysis will form a mathematical extension of the autocorrelation 
function, through Fourier transformation, to yield the power spectral 
density function. 
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AIRCRAFT DESCRIPTION AND FLIGHT TEST PROCEDURES 



The basic measurement used to determine atmospheric turbulence 
In flight Is the angle of attack of the aircraft as measured by flow 
vanes. The flow angles measured are corrected for the effects of 
airplane motion, Chilton (3). In this particular analysis as applied 
to the power spectrum method, a continuous time history of the vertical 
component of the trie gust velocity Is to be Initially determined. 
Although the lateral and longitudinal components of turbulence could 
have been measured as well, the vertical component has the most slgnl* 
f leant effect upon the aircraft's pitch dynamics, and In the Interest 
of simplification only this vertical portion of the turbulence will be 
considered. 

The vertical component of the gust velocity, from Ref. 1, Is 
given by: 

The equation above Is based upon the following assumptions as outlined 
In Ref. 1: 

1. All disturbances are small enough to allow the use 
of the angle in radians In place of the sine of the angle. 

2. Boom flexures are negligible. 

3. The effects of variation In upwash on vane Indica- 
tions are negligible or allowed for by calibration. 

In conducting the flight analysis, the average values of the 
measured parameters (denoted by a bar) are Included to establish the 
flight measurements as increments from the mean values for the measured 
flight Interval. Therefore, the time history of the vertical compo- 
nent of the gust velocity, from Ref. 1, is: 

t- 

A reading Interval of 0. 1 seconds was established as being an 
optimum value for this analysis after considering the aircraft speed at 
which the data were obtained (300 knots TAS) and the frequencies 
evaluated . 
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All turbulence date used £or this analysis were obtained by 
the author during the month of June 1967 while participating in an 
industrial tour at the Flight Test Branch of HASA's Ames Research 
Center at Moffett Field, California. The flight test data used for 
this thesis are only a portion of the total data obtained by NASA as 
part of a larger program to study the optimum dynamic characteristics 
and response to turbulence of a variable stability F-IOOC airplane. 

The aircraft used for this turbulence study and its principal 
dimensions are shown in Fig. 4, from Ref. 4. It is a North American 
F-IOOC which is assigned to the National Aeronautics and Space Adminis- 
tration and designated as NASA aircraft' number 703. The F-lOO is a 
single place, low wing aircraft, powered by a J-57 aeries turbojet 
engine with afterburner. The wing and horizontal stabilizer are swept 
back at an angle of 45 degrees. Low speed flight and stall character- 
istics are improved by the use of longitudinal leading edge slats on 
the standard F-lOO airplane, and good longitudinal control is provided 
by an all-moving "slab" tall. The controls are completely power 
actuated. Longitudinal stick forces are provided by a simple spring 
bungee whose leverage on the stick is constant at all speed ranges 
except between a Mach number of 0.8 and 0.94 where the leverage is 
altered by a Mach number sensor. The F-lOO airplane has a 400 cycle, 
115 volt A.C. electrical system and a separate 28 volt D.C. system. 

The parameters measured in flight are given by Equation (5), 
and the type and location of the flight instruments are shown in 
Table 1. A continuous record of aircraft response was necessary in 
order to provide a complete and effective record of the turbulence en- 
countered, therefore an oscillograph. Instead of a photopanel, was 
required. All data were recorded on a 26 channel photographically 
recording oscillograph and developed on Ektachrome Aero film for 
analysis. As can be seen in Fig. 5, the accelerometers were positioned 
as close to the aircraft's center of gravity as possible. The flow 
vane used to measure." angle of attack was located on the top of the 
vertical stabilizer as shown. A standard pitot static head was 
utilized for airspeed measurement. It is located ahead of the aircraft 
on the tip of the F-lOO's standard nose boom. 
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AIRPLANE GENERAL INFORMATION LO. 1F-100C-2-5 

AIrpian# Dimmsioni 





1-12 




- VO di (5) 



PARAMETER 

MEASURED 



TYPE OF raSTRlTrlENT 



LOCATION 
OF IMSmUIffiNT 



INSTRUMENT 
LffilTS 



Airspeed 

\/ 



Standard F-lOO pitot- 
static instrument with 
Bourns pressure cell 



Nose boom 



0-^00 kts 



Angle of 
Attack 

OCy 



IlOO cycle A.C. NASA 
instrument 



Top of vertical 
stabilizer, 
Station 52^ 



t 10° 



Pitch Angle 

a 



Mirneapolis Honeywell 
attitude gyro with 
potentiometer pickoff 



Left link 
compartment 
Station 1^2 






Vertical 

Acceleration 

a. 



Statham linear strain 
gage accelerometer. 
Model CS-8-350 



■flJheel well on 
aircraft's 
longitudinal 
centerline 
Station 315 



t. lig's 



Pitch Rate 

O 

e 



Spring restrained rate 
gjrro with a Micro syn 
signal generator 



Left link 
compartment 
Station 114.3 



Jt90 / sec 



Longitudinal 
Stick Position 



Potentiometer with itOO 
cycle A.G, output 



Stick linkage 
on centerline 



:^-19°stick 

deflection 



Horizontal 

Stabilizer 

Position 



Linear ootentiometer 
pickoff using Mar kite 
infinite resolution oot. 



Station 



^ 25 °, - 5 ° 

stabilizer 

position 



Altitude 

h 



Statham pressure cell 
with strain gage pickoff 



Nose compartment 
Station 55 



0-50,000 ft 



TABLE I: AIRCRAFT INSTRUMENTATION: NASA AIRCRAFT NUMBER 703 
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This particular F-IOOC had been modified by NASA for use as a 

variable stability airplane in order to study the effects of the 

interaction between the turbulence and the longitudinal stability 

characteristics « A planned variation in C and C was developed 

m 

by NASA engineers, and these two parameters were ^ changed to pro- 
vide varying degrees of pitch response. 

From basic theory, the values of C primarily determine the 
natural frequency of the short period mode and are a major factor in 
determining the response of the airframe to elevator motions and gusts. 

The stability derivative, , is usually referred to as the 
pitch damping derivative and is ^ critical in longitudinal dynamics 
because it establishes a major portion of the damping of the short 
period mode. This damping effect comes mostly from the horizontal 

stabilizer. The appropriate C and C gain settings were obtained 

mo6 m 

by means of an analog computer matching ^ technique to produce the 

desired combinations of the short period frequency and damping ratio. 

For this analysis, by varying C and C , a resultant value 

of equal to -14.47 and M equal to .3068 ^ was programmed into 

^ q 

the aircraft to produce a short period frequency of .6 cycles per 
second and a damping ratio equal to 0.10. 

The simplified schematics in Figs. 5 and 6 show how the variable 
stability mode of this NASA F-IOOC is obtained. Three modes of control 
may be manually selected by the pilot: a. Basic Mode, b. **Fly By 

Wire** Mode and c. Variable Stability Mode (includes **Fly By Wire** 
mode) . 

In the variable stability mode as shown in Fig. 5, a signal 
from both the angle of attack vane and the pitch rate gyro are received 
at the summing amplifier. The output from the summing amplifier is 
summed with the pilot's input and fed into a servo which controls the 
horizontal stabilizer* s hydraulic actuating valve to provide the 
desired response characteristics. 

A detailed schematic of the pilot *s input as a closed loop 
system is shown in Fig. 6. A potentiometer is used to detect control 
stick inputs with a sensitivity such that the basic aircraft stick-to- 
stabilizer gearing is preserved. A torque tube between the stick and 
the hydraulic valve is normally mated firmly in the Basic Aircraft Mode. 
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FIG. 6 SIMPLIFIED SCHEMATIC OF CLOSED LOOP CONTROL 
SYSTEM IN NASA F-IOOC AIRCRAFT NO. 703 





In the "Fly By Wire" or Variable Stability Modes, the torque tube is 
mechanically adjusted by the system to have a greater increment o£ 
play. This adjustment prevents any feedback from the stabilizer to 
the stick; and allows additional comsiands to the stabilizer in the 
Variable Stability mode. 

C and C are altered by the variable stability system in 
accordance with the following equations: 



0/r?, 






= c. 



m. 



^oc 6a S/C 









(7) 



oc. 



C/77^ eff Sas/c 







( 8 ) 



In each example, the basic values of the F«100's stability 
derivatives are modified by the second term in the equation whose 
value is altered by the servo system. 

Since the longitudinal response of the F-lOO airplane to tur- 
bulence was of prime importance, all flight tests were conducted with 
no longitudinal control inputs, however the wings were kept level, 
while data were being taken, by applying minor lateral stick inputs. 

As a normal procedure, the turbulence was entered from a stabilized 
flight condition at an airspeed of 300 knots true, which is a realistic 
cruise airspeed for the F-lOO airplane at the test altitudes selected 
(13,000 to 15,000 feet). Although no longitudinal control inputs 
were desired, both stabilizer position and control stick deflection 
were checked continually when analyzing the flight test data to insure 
the validity of the measured aircraft response. 

Most of the turbulence encountered was of a mountain wave or 
clear air turbulence type. The leeward side of Mount Whitney and the 
neighboring peaks in the Sierra Nevada range were used as a flight 
test area. 
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BASIC COMPUTATION PROCEDURE 



As in any statistical approach to analyzing data and predicting 
trends, the more information obtained, the more valid' the study. This 
concept was tempered by the selection of that flight test data which 
were obtained with the F-IOOC programmed for the most realistic short 
period frequency. Accordfpgly, after screening the available test data 
NASA Flight 26, Run 8 was selected as being a valid and optimum flight 
record for the power spectrum analysis. On this flight, the F-IOOC 
was programmed for a short period frequency of .6 cycles per second. 

It is reiterated that the intent of this thesis is to demonstrate the 
application of the power spectrum analysis to determine the dynamic 
response of an aircraft to turbulence rather than to establish any new 
information on the F-IOOC airplane. 

The data for this flight were read every 0.1 seconds, and the 
basic parameters as outlined in Equation (5), were recorded on stan- 
dard IBM eighty column digital computer cards. Due to the large quanti 
ties of data Involved (six inflight parameters read at 920 intervals 
each) and the numerous computations required by the power spectral 
method, a digital computer is necessary to effectively process the 
data and perform the analysis efficiently. The average values of the 
measured flight parameters, 'denoted by a bar in Equation (6), were 
initially determined by means of a simple digital operation, and w 

8 

itself was computed. The power spectral analysis of the measured 
flight data was performed by a digital program utilizing an IBM 
7094/7040 direct coupled system computer for the calculations and 
graphing. 

Both the autocorrelation function and the power spectral den- 
sity of the input, w^ (referred to as series 2), and the output, & 
(referred to as series 1), were determined and graphed by the computer 
program. The desired transfer function, was then computed in 

both phase and amplitude as a function of gust frequency. The transfer 

function was based upon a cross correlation and cross spectral density 

« 

function which combined the average values of Q and w over a 

s 

common frequency interval. From the computer analysis, the phase and 



amplitude of 



were graphed on a cross spectral density plot 
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versus frequency. To check the accuracy of the power spectrum method, 
the theoretical solution for the phase and amplitude of r-— was 
obtained by solving the longitudinal dynamic equations of motion and 
the pitching moment equation about the lateral axis. The phase and 
amplitude of the transfer function obtained theoretically were also 
plotted versus frequency, and the resulting curves from both methods 
were compared. 
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DISCUSSION OF THE POWER SPECTRUM METHOD 



A detailed analysis of the steps involved in the digital com- 
puter program is presented in Appendix A; however, a proper apprecia- 
tion of the digital computation can only be realized if the mathe- 
matical theory behind the power spectrum method is outlined in detail. 

As discussed earlier, the vertical component of the atmospheric 

turbulence, w , was established as a stationary, ergodic form of 
S 

random data. An additional stipulation is that the turbulence should 
represent a Gaussian process. For the turbulence to be Gaussian, the 
fluctuations of the turbulence must have a normal probability distri- 
bution. References 1 and 5 stipulate that if one component of the 
turbulence, w , for instance, has a normal distribution, then normality 
of the turbulence can be assumed. Figure 7, which is reproduced from 
Ref. 1, shows the cumulative frequency distribution of w^ as determined 
by NASA from a thunderstorm traverse. The plot is scaled in such a 
manner that a normal distribution of w will result in a straight line. 

s 

The data were obtained from a 200 second record, read at 0.1 second 

intervals. A fitted normal distribution is shown in comparison to 

the data points plotted, and it may be noted that there is a good 

correspondence between the fitted normal distribution and the traverse 

data points, except at large values of w , therefore, from Ref. 1, 

S 

the assumption that there is a Gaussian distribution of turbulence 
appears to be reasonable. It is reiterated that the properties of 
statlonarity, homogeneity, and isotropy, in addition to the Gaussian 
assumptions, are Implidd. These properties specify, then, an invariance 
in the statistical characteristics of turbulence with respect to the 
position along the time history of the sampling interval and the 
direction through the turbulence. 

From Ref. 6, Taylor's hypothesis that the variation in gust 
intensity existing along the flight path at any given Instant will 
remain substantially the same until the aircraft has traversed the 
given body of air, is utilized. The relative length of the body of 
air along the flight path, therefore is large compared with the dis- 
tance traversed in the reaction time of the aircraft. 
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Vc/'i- /'CQ./ Co/77pone/-)i~ o-f Qus-f" \/<s/oc/i~y^ ff/sec. 



Figure 7. - Probability of equaling or exceeding 

given values of vertical gust velocity 
for thunderstorm (traverse 4). 
(Reproduced from Ref. 1) 



36 



with the mathematical assumptions established, a detailed 
analysis of the mathematical steps behind the power spectrum analysis 
can now be presented in conjunction with appropriate sequential 
reference to the computer program which is detailed in Appendix A. 

A block diagram (Fig. 17) showing the steps involved in the 
computer analysis and the appropriate mathematical expression for 
each step is presented in Appendix A in order to provide a complete 
picture of the entire digital program involved in the power spectrum 
analysis. 

The magnitude of O and w , as computed from Equation (6), 
is tabulated at 0.1 second intervals on the computer printout and 
referred to as series 1 and series 2 respectively. The values of the 

data are read from left to right in each row in chronological order. 

« 

The autocorrelation function of O and the autocorrelation function 
of w are computed separately by taking the product of S or v at 
time t and at time t 7 , and averaging these values over a total 
observation time, T, The resulting average product will approach an 
aacact.aistO0orr elation function (also referred to as the autocovariance 
function) as T approaches infinity* 

In equation form the autocorrelation function is expressed as; 



r 

^ f' X (t) • X (t y) c/t 

- 7 — — OQ 



(9) 



where X(i) represents the value of either & or w^ at time t, and 

/(ti- y) is the value of ^ or w after a time shift whose magnitude 

S 

is determined by the value of 7 . 

s 

The autocorrelation functions of & and w^ are each computed 
for progressive values of y , that vary in 0.1 second increments, 
from 0.1 seconds to 9.99 seconds. The autocorrelation function for 

i 

0 is tabulated on page 77 of Appendix A and plotted versus time 
on page 78 also* Similarly tabulated on page 82 

and plotted versus time on page 83 in Appendix A* 

The quantity , as defined above, is a real valued function 

with a maximum value at 7 • 0* ^Iso and 

for all 7 * The mean value of x(t) can be determined from » 

i.e., it is equal to the positive square root of the autocorrelation 
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function, as the time displacement increases toward infinity (Ref. 2). 
In addition, the mean square value is equal to the autocorrelation 
function at zero time displacement, (7 “0); i.e., (^o) 

As an example of the usefulness of the autocorrelation function, 
a simple autocorrelogram (plot of an autocorrelation function) is 




(Reproduced from Ref. 2) 

Although phase information is not available with this plot, it 
can be noted that the autocorrelation function for the sine wave will 
persist over all time displacements since it represents deterministic 
data. Conversely, an autocorrelogram for narrow band random noise, 
(non>deterministic data), will diminish to zero for large time dis- 
placements, assuming a mean value of zero. This is illustrated in 
Fig. 9 




o 

Figure 9. - Autocorrelogram for Random Noise 
(Reproduced from Ref. 2) 

A combined plot of the autocorrelation function of the sine 
wave superimposed with the autocorrelation function of the random 
noise will appear as in Figure 10. 
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Figure 10. - Combined Autocorrelograms for 
a Sine Wave and Random Noise 
(Reproduced from Ref. 2) 



The deterministic data (sine wave) can be detected from the 
random background of the wide band noise. This is the principle 
value of the autocorrelation function in that the random and non- 
random data may be separated into their respective categories through 
analysis of their combined autocorrelograms. 

o 

The autocorrelation function of ^ , series 1, is noted to 

diminish with increasing values of time, due to the inherent longi* 

tudinal damping characteristics of the F-lOO airplane. This is a 

typical response of a damped second order system to a random input. 

Refer to pages 77 and 78 of Appendix A. 

It is also to be noted that the autocorrelation function of 

w , series 2, is also a maximum at 7 » 0 and decreases to a minimum 

8 

value for increasing values of 7 (pages 82 and 83 , Appendix A). 
This is the expected characteristic of the autocorrelation function 
for random data. 

After establishing each autocorrelation function, and 

f the power spectral estimates of both ^ and w are computed 

V 8 

separately and tabulated on pages 79 and 84 of Appendix A. The 
power spectral estimate is simply a Fourier transform of the auto- 
correlation function and describes the general frequency composition 
of the data in terms of the spectral density of its mean square value, 
(Refs. 2 and 7). The time tbean square value provides a 
measure of the disturbance energy per unit time and thus has been 
referred to as power in the study of turbulence and has been termed 
the energy or power spectrum. The power spectral density function 
specifies the contribution of the harmonic components to the mean square 
value of the random data over the specified frequency range. 
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The average squared value of the random data will approach an 
exact mean square value as the observation time, T, approaches infinity 
for a frequency range between f and ^ f (Ref. 2). 

In equation form: 






( 10 ) 



T'—^oo 

where x(t,f, A. f) Is that portion of x(t) in the frequency range 
bounded by f and f + ^ f * In this analysis, of course, x(t) could 
represent either ^ or v , For small values of f, the power 

s 

spectral density may be computed to be: 



( 11 ) 



In precise form; 



^ '' ( 12 ) 

O ' 

where G (f) represents the power spectral density function of either 
• ^ 

^ or w . The power spectral density is always a real valued non- 
S 

negative function. 

As noted earlier, the power spectral density is a Fourier trans- 
form of the autocorrelation function and for stationary data may be 
expressed as: 

C30 



<?o 



4/k^(y)coS ztt/ 7cf7 



(13) 



Thus, from Ref. 2, the mean value of x(t) is given by ^ 
mean square value is given by ^ 

The Siean square value, therefore, is equal to the area under the plot 
of the power spectral density function versus frequency. 

From Ref. 1, typical power spectra for the vertical component 
of turbulence in clear air, cumulus clouds and thunderstorms are 

» 

plotted in Fig. 11. The power spectral density is in units of (ft/sec)' 
per radian/ft, and it is plotted versus the reduced frequency in 
radians/ft, with an additional scale of wavelength in feet, also along 
the abcissa. Since the gust disturbances are essentially space 
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Power speciral density, ({>(11), (ft/sec) Vradlans/ft 



100,000 




Reduced frequency, , radians/ft 



Figure 11. Typical power apectra of vertical component of 

turbulence measured in clear air, cumulus clouds 
and thunderstorm. 



41 



disturbances rather than disturbances In time, a reduced frequency, 

.y2. > in radlans/ft. Is normally used In the power spectral analysis 

of turbulence. The relative heights of the curves In Fig. 11 Indi- 
cate the variation of the Intensity of the turbulence encountered 
under the different weather conditions. The similarity of the slopes 
of the power spectra Is apparent, and their magnitude Is equal to 

, where C Is a constant. The spectra show a characteristic 
rapid decrease In power with Increasing frequency. 

The root mean square value of the gust velocity Is a measure 
of the Intensity of the turbulence and Is equal to the area under the 
portion of the spectrum evaluated. From Fig. 11, it may be noted 
that the root mean square value of the turbulence varies In stagnitude 
from 4.48 ft/sec in clear air to a value of 13.38 ft/sec In a severe 
storm. The above example is presented to describe the type of Infor- 
mation that is nonully available from the power spectral density 
plot of a given parameter. 

O 

The power spectral density of & and w are plotted on pages 

8 

80 and 85 of Appendix A. As expected, the power spectral density 

of w (series 2) diminishes with Increasing frequency. The power 

8 • 

spectral density plot of Q (scries 1) also decreases with increasing 
frequency; however, a significant mode of response at the programmed 
short period frequency of .6 cycles per second Is apparent. 

e 

The general dependence of one set of data (such as d ) upon 
another (w^) auy be described In more detail using the cross correla- 
tion function. From Ref. 2, two time histories of data are Illustrated 
in Fig. 12. 
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Like the autocorrelation function, the average product of x(t) 
at time t and y(t) at time t + 7 is taken over the observation time 

a 

t. For this analysis x(t) is comparable to 0 » and y(t) is comparable 

to w or vice versa. As T becomes infinite, the resulting average 
S 

product approaches an exact cross correlation function. As illustra^ 
ted in Ref. 2: 



y^)cyCy)= 



(14) 



R ( ) is always a real valued function which may be either posi* 

xy 

tive or negative; however, unlike the autocorrelation for random data, 

'^yy('y) may not necessarily have a maximtim at 7 ■■ 0 or be an 

even function only. R ( C7 ) displays S3nnmetry about the ordinate 

xy 

when X and y are interchanged, so that R (7 )"R (-ZT). 

xy xy 

The absolute value of the cvess correlation function is bounded 



by: 



/ f y)/^ - (o) ' ('o) 



/ 



//?^y (y)/ ^ i H (O) ^ (O)/ 



where. R is comparable to the autocorrelation function of w , R is 

X « 8 y 

comparable to the autocorrelation function of 0 , and R^^ is the 

cross correlation function based on 4 and w . The cross correlation 

.8 

function (labeled as cross covariance) for A and w is tabulated on 

8 

page 86 of Appendix A, for both a positive 7* add a negative 7 , 
and plotted versus time on pages 87 and 88. 

The Fourier transform of the cross correlation is taken to 
determine the cross spectral density function which is comparable to 
the power spectral density function for a single variable. Since the 
cross correlation is not an even function, the cross spectral density 
function is a complex number composed of the co-spectral density func- 
tion (real part) and the quadrature spectral density function (imagin- 
ary part). Therefore: 



(f) = Cxy (f) - -^'Qxy ( 15 ) 
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where !• the co-spectrel density function and is the quadra- 
ture spectral density function. 

The co-apectral density function is the average product of 
x(t) and y(t) within a narrow frequency interval between f and f 
^ f divided by the frequency interval. 

The quadrature spectral density function is sinilar to the co- 
spectral density function, except that either x(t) or y(t) is shifted 
in tine to produce a 90 degree phase shift at the frequency f. 

The co-spectral density function is expressed as: 



/ 

o 

and the quadrature spectral density function is: 



(16) 






(17) 



%rhere x(t,f,^ f) and y(t,f, ^ f) are filtered portions of x(t) and 
y(t), and y*(t,f, A f) signifies a 90 degree phase shift from y(t,f,^f). 
The co-spectral density function and the quadrature spectral density 
function for Q and w^ are tabulated on page 89 of Appendix A. 

The cross spectral density function nay also be expressed in 
complex polar form as: 



The magnitude/G (f)/ of the cross spectral density function is 
xy 



equal to: 



/G,y (/)! = dc^y(/> ^ Q\y(d 



(19) 



and the phase of the cross spectral density function Q (f), is 

xy 

expressed as: ru v 

( 20 ) 
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The magnitude and phase of the cross spectral density function 

are tabulated on page 89 of Appendix A and plotted on pages 90 and 

91. Again a significant peak in the amplitude of the cross spectral 

density function is noted in the region of the programmed short period 

frequency of ,6 cycles per second. 

The amplitude of the transfer function, between the input, 

x(t), and the output, y(t), is obtained from the amplitude of the 

cross spectral density function, G (f), and the amplitude of the 

xy 

power spectral density of the input, G^(f). The amplitude of the 
transfer function, input x(t) to output y(t), is expressed as 
and is equal to: 







In this analysis x(t) corresponds to w 
so that: 



g 



and y(t) 




ft 

is likewise 0 



( 21 ) 



( 22 ) 



The phase of the transfer function is the same as the phase 
of the cross spectral density function as outlined above. 

Q 

The amplitude and phase of versus frequency are tabulated 

on page 92 of Appendix A, where the p^ase of the transfer function is 
expressed as portions of a circle. The appropriate translation to 
angular magnitude; in degrees, is noted on the computer program beside 
the printed figures. 

Q 

Two separate curves showing the phase and amplitude of 
versus frequency were plotted from the results of the power spectral 

analysis and are included as Figs. 13 and 14. It is noted that the 

& ^ 

peak amplitude of , and also a phase shift of , occurs at the 

programmed short period frequency of .6 cycles per second. Although 
scatter of the data points is apparent at frequency values greater 
than the short period frequency, the amplitude peak of the transfer 
function and a distinctive phase shift may be readily determined from 
the plotted curves. 
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THEORETICAL SOLUTION FOR THE PHASE AND AMPLITUDE OF THE TRANSFER 

FUNCTION 



In order to determine the accuracy of the phase and amplitude 



6 



of ^ , as computed by the power spectrum method, the appropriate 
longitudinal equations of motion for the F-IOOC airplane were solved 



on a theoretical basis. 



Using the programmed values of C^^ and 



0^, as well as the standard stability derivatives for the F-IOOC; 
the phase and amplitude of were computed over the same frequency 
range as that used for the power spectrum analysis. 

The following basic assumptions were made in determining the 
equations of motion and solving for — : 

1. The equations of motion were written with respect to the 
stability axes. 

2. The aircraft’s elevator was considered to be locked, so 
that there were only three degrees of freedom. The equa- 
tions were written therefore for linear motion along the 
X and Z axes and for pitching moments about the Y axis. 

3. The aircraft was considered to have a linear response to 
the turbulence encountered. 

4. No lateral/longitudinal coupling occurred. 

5. The airspeed was stabilized at the beginning of the run 
and the engine thrust was constant throughout the run. 

6. Quasi-steady flow was assumed. 

7. Only the vertical component of turbulence, w , contributed 
to the aircraft’s normal acceleration. 

The phase and amplitude of the transfer function between w 
• 8 
and S were initially obtained using a simplified solution in which 

the phugoid motion was considered to have a negligible effect upon the 
short period oscillations due to the turbulence. See Ref. 9. This 
simplified method is designated as method A in this section. To check 
the theoretical results, a second solution involving all of the essen- 
tial terms included in the stability equations is included. This exact 
method is designated as method B. Both data sheets for methods A and 
B are included in this analysis together with the appropriate values of 
the stability derivatives and the important terms used in each solution. 
The results of both theoretical solutions were in close agreement. 



The phase and amplitude of obtained by methods A and B are 

plotted in Figs. 15 and . 16 . It can be noted that the peak amplitude 
and the phase shift of — occurred at the programmed short period fre- 
quency of .6 cycles per second. 

From Section 13 of Ref. 9, the longitudinal equations of motion 
are developed with respect to the stability axes: 



Xax/s: d r ge cos = Tu (coS f)u + cosf 

V- i- X<^q i- X^Wf Sr -h XifrSr X^^Sa 

XrOK'^ '- S/n ^ 

^tchuo. g -- ^ 



(23) 



(24) 



(25)' 



These equations may be further simplified since the angle Ot, is 
assumed to be small so that sin 0£> » 0 and cos “ !• In addition, 
since there is no elevator deflection, flap deflection, or speed 
brake deflection during the run, Si£ ^ S p- and Sg all equal zero. 

Also under the Initial conditions, the thrust is assumed to be con- 
stant; therefore, T^ is equal to zero, and the variation of the thrust 
with RPM, T . is also zero. X , the variation of the force along 

the X axis, and Z , the variation of the force along the Z axis, due to 
<1 

the angular pitching velocity, are negligible. In addition X and Z , 

w w 

the variation in the forces along the X and Z axes due to the time 

rate of change of the vertical velocity, are negligible. For this 

analysis, the vertical velocity is considered to be composed of two 

components, w and v , where w is the steady state component of the 

S 

vertical velocity and w^ is the fluctuation of the vertical velocity 
due to turbulence. 

Eliminating the terms as noted above, separating the vertical 
velocity into the steady state and turbulent components, and intro- 
ducing La-Placiar notation, the equations of motion become: 
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5 U^<je = x^u ^ X^u + Xi^Wg ( 26 ) 

^ -c/ose = f tv ^ 

s^e^MuU-^ /^^se ^ 

Re'-«rranging the equations and solving for the w^ terms: 

(s-xu)^ -x^y -f- c^e - < 26 ) 

-Zi,u -sOoO = £tv -< 27 ) 

-Mc(U -(xi^S^M^)i^-^S(s-Xl 2 ^ 0 =/^^yg ( 28 ) 



Using Cramer's rule, the above determinant Is modified to solve 
for so that the transfer function Is equal to: 





(s-Xu) ~)(w 


Xi^ 




-Zu (5- ?w) 




e _ 


,-Alu 






(It-Xc/) “Xiv 


? 






-J0> 




- Alt< S i- XI J) 


S(S-M^ 



( 29 a) 



Expanding the determinant: 

-AAufJ]-hM^f(3- X^)(S-?w) - Bu Xvj] 
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where ^ denotes the original determinant in the denominator o£ 
Equation (29a). Equation (29b) can be further simplified to: 

-f Xw Mu 

^ ^ (M^ -Xu Alt;,) - fXu ^ ^^v)Js (29c) 

Multiplying the numerator of the ratio of the determinants by s, an 
expression for is obtained, so that: 

^ (^W -XuM^) - Mw (Xu / ^1^)1 5^ ( 30 *) 

Method A: The Simplified Solution of the Longitudinal Equations 

of Motion 

For the simplified solution, the original determinant of the 
denominator is considered to be equivalent to the product of the 
quadratic expressions for the short period and phugoid motions, in 
Lapliicdixi : form, so that: 



^ J 1 / 

where and are the damping ratios of the phugoid and short 

period motions respectively. And (*->pp and uj^p are equal to 2 7T 
times the phugoid and short period motion frequencies. The contribu- 
tion of Ldpp is considered to be of negligible significance in the air- 
craft’s dynamic response to turbulence so that the terms containing 
can be eliminated. Therefore, the denominator can be approxi- 
mated by: 



‘■Xs^ ^ 2 OJ^p 3 + CU^) 

2 

Canceling s in both numerator and denominator and replacing the 
operator with jw- , the simplified expression for becomes: 

4 — ^ ^XwMlL ^ — ()( u 

Oj"^) ^ (32) 
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The linear velocity along the X axis is considered to be constant so 
that the perturbations Z^, and may be neglected. Equation (31) 
therefore reduces to: 



^9 (oJsp - oJ^) -^'2 (^sp cJ 

The amplitude of the transfer function (method A) is; 



0 

h/^ 



and the phase 



^(ujJ-p - 4 cOsp ' 

of the transfer function is: 



(34) 




90 °— -tan"* Z uJsp cj 

(Co ^p — to ^3 



(35) 



The expressions for the stability derivatives are given by 
Ref. 9 as: 



M 



^ -_£SV /c, ^ Q,) 

2/T7 ^ ^ 



A1 






w 



Also from Ref. 10: 



^4/ 









-Zoc - Mi 



(36) 

(37) 

(38) 



(39a) 






where: 

z 



(40a) 






- a c.;) <”> 
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and therefore: 



~ ^sp ^^pJ 



(40b) 



A/^ -- 




(39b) 



The computed values of the stability derivatives noted above, 
as well as the quantities used to compute them, are listed in Table 
11. Using the values of the stability derivatives as computed in 
Table 11, the phase and amplitude of are computed for a range of 



frequencies from .05 cycles per second to 2.00 cycles per second. The 
values of the phase and amplitude of , as computed by the simpli- 

fied solution (method A), are listed in Table 111. The phase and 
amplitude of , as computed by the simplified solution, are plotted 

versus frequency in Figs. 15 and 16 . 



Method B: The Exact Solution of the Longitudinal Equations 

of Motion. 

The initial computational steps of the simplified solution, 
method A, as well as the same given stability equations and initial 
assumptions are implied for the exact solution as well; however, for 

the exact solution the entire determinant for w will be retained in 

° g 

the denominator of the term throughout the analysis. Starting 
with Equation (30a) , ^ is equal to: 

where is equal to the expansion of the determinant formed by 

the longitudinal stability equations, solved for w . 
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TABLE II: FLIGHT DATA AND COtPUTED STABILITY DERIVATIVES 
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TABLE III: METHOD A. SIMPLIFIED SOLUTION FOR THE PHASE AND AMPLITUDE OF 
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TABLE III; (Continued) 
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^ = Sfs-Xu)fs- ^)(s- - S(Jo Cs-XuX'^i^S^ ^i^k9iafx}^S^/lX) 

- J- A^v' ^cj f J- A'/^) -^gA1u{s-2r^)-S Uo X^ XI u <^3b) 

A = s"^— (Xl( ZiV a a Uo M^^)S'^-^ ^X'u A ^(v) A Xu ? ^~Xu>iu. 

- Uo [Uc (XuX^kJ “/v Mu)-^XI^(Xi^^a ~Xq Xtv) 



^ A ^/^^ 6 /A/iv — A^a i?»v) 



(43c) 



For simplification, the real and imaginary parts of the numera- 
tor and denominator are considered separately and designated as: 



*4) 



The real portion of the denominator 
The Imaginary portion of the denominator 
The real portion of the numerator 



N 



The Imaginary portion of the numerator 



expressed as: 



o 

& 

Therefore, the amplitude of the transfer function - 7 — 



^9 1/ A 



and the phase of the transfer function is: 



is 



(44) 



^ ^ = —Xon ' X'9 — fa/'/ ' 

/?v ^£> 

where: 

Aw o, ~X^Mu 

^ 10 ^ A (46a) 

7^ - W (47) 

A/, • ^ f> (47a) 

J ^ - y '{ ty/ A 7 ^ ■/ Uc A! Vs/) A Uo (Xu f^v/ — Xh/ AJu) 

A /'^l^ (Xi'v' f ~ X.i jf 1^) /■ y^Alu f u XI A')J' 

_7^ c.(of/: / '■- ) (^ 8 a) 
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(49) 



- Uo (M^ - Xc( M w)] - Mu 






Cu"^ - CaJ^C 



(49a) 



The stability derivatives as listed for method A, and their 
values as tabulated in Table II, apply to method B as well. In addi- 
tion, the following derivatives, computed in accordance with Ref. 9, 
are used to solve for the coefficients A through F, of Equations (46a) 
through (49a) : 



-- C 



rn 



£> 



(scH, 



'^L 

^C( - ^ 



(52) 



(zTy^O^^<^ < 53 ) 

The computed values of these stability derivatives, as well as 
the quantities used to compute them, are listed in Table II. The values 
of A, B, C, D, E, and F are listed in Table IV. The phase and amplitude 



of 






are computed and listed in Table IV for a frequency range from 



.05 cycles per second to 2.00 cycles per second. The phase and ampli- 

tude of -X— , as computed by the exact solution, arc also plotted versus 

h/o 

frequency, together with the values obtained from method A, on Figs. 15 

and 16. The two different theoretical methods produced almost identl- 

A 

cal curves. For both methods the peak amplitude of 77- , as well as the 

h/g 

phase shift, occurred at the programmed short period frequency of 0.6 
cycles per second « 
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TABLE IV: METHOD B. EXACT SOLUTION FOR THE PHASE AND AMPLITUDE OF 






■V. 


K 

\ 


(S 

VS 

\ 

\ 


NS 

K 

\ 
1 * 


CK 

V 
\ ♦ 


Os 

53 

N 

Ni 

in: 


N 

R 

iS 


|N 


>0 


N4 


1 


1 

•'i 

1^ 


§; 

1^ 


?? 

1R 


Os 

M 

f>l 

t N 


1 

N| 

>: 

|N 


,>i 


o; 

IN 


Ck 

‘h 

^'^4 


1 

1^ 


1 

,« 


N 

1^ 


5 

,R 




•N 

N 

S 




N 

a 


0 


Os 

§ 

l'? 


>s 

,1 


K 

N 


N 

'?: 

K 

N 

IS 


«K 
\ < 


Os 

1 n‘ 


s 

K 


1 

I^ 


Ns 

1^ 


s 

1 K 


On 

Ns 

5: 

Ns 

«< 

1 


SI 


N 

V 

N 


15 

<k» 

V 

I'N 


$? 

N 

O 

N 

l'^ 


K 

V* 

-s: 

0^ 

IN 


R 

-N 


1 

5 


fNj 

g 


'5 


N 

K 

,K 


1 

1 


1^ 




V 


o 

N 


CK 


1 

iS 


C5 

<K 


K. 

§ 

In 


i 

6^ 

1^ 


K 

r^ 

fS 

t5 

K 


Nft 


\ 

1^ 


1 

1^ 


N 

si 


s 

Ss 

K 


V^ 

0 

v» 

Vi 
1 . 


s 

N 

v> 


1 


N 

N 

N 


\fi 

It: 

N 


Os 

N 


)Q 

SJ 

K 

Si 


t 

B 


Vi 

1 


1 


1 

\ 


1 

1 




V' 

.3 


O 


s 


s 

3? 

iS 




N 

K 

a 

On 

IN 


if? 


K 

V 

? 

K 


1 


Nv 

K 

•o 


N 


R 

K 


K 

tN 


<3 

Ns 

\ 

VI 

v> 
1 • 


Vi 

K 


5? 

Os 

Os 


•n 

1 

§ 


1 

1 


fN 

Ci 

N 


s 

Ns 


Sj 

i^ 

1 


1 

1 


1 


«6 

1 

^ 


N 

Os 

N 




V 

N 

3 




N, 


CK 

N 

"s 

•N, 

v» 

s 


§ 

N* 

N 






r4 

vS 


1 

5$ 


K 


Os 

N* 


§J 

Os 

s 


Os 

CS 


V 

fs 

Vi 

5 

-A 


Vi 

s 

Vi 

K 

c? 


Ns 

S5 

<K 

V 

N 

N 


>« 

1 


1 

1 


f 


K 

Ns 

n; 

4 


<3 

On 

Os 

v> 


o: 

On 


§ 


tN 

M 

N 

N 

1 


M 








S 

O 
1 • 


's. 
1 • 


On 
1 ' 


oo 

Oo 

> 

f • 


VJ 

\ ' 


00 

\ 

1 ■ 


1 
N 
1 • 


K 

-5 
>0 
s 
1 ' 


V> 

N 

■S. 
1 ' 


o 

•r> 

•S. 

0 
N 

1 ' 


1 • 


N 

•'< 


? 

tS. 


'IK 

s? 

,N 


;s 

Vi 

iN 


»s4 

!^ 

Oo 

|N 


fc- 

On 

oo 

Vi 

1^ 


1 

Ci 

Ns 

VI 

y> 


V 

On 

VA 


Os 

V 

On 

1^ 


r>l 

Vi 

Os 

N 

1*^ 


Qo' 

On 

Ns 

K 

N 

I'n 


N 

1 

;s 






0 


N 

§ 


(K 

O 


0 


<0 

§ 


V 

V5 


Ci 


N^ 

Ci 

Ns 

0 


On 


V® 

S 

S, 

Vi 

Vi 


K 

Ns 

N 

Vi 

Vi 


*Vi 

N 

Vi 

Vi 


'>3 

Ns 

Vi 

Vi 


K 

M 

Vi 

Vi 


VS 

NJ 

V 

Vi 

Vi 


Os 

V 

Vi 

Vi 


1 

Vi 


Vi 

Vi 


VO 

Vi 

Vi 


\ 

Vi 

N 


N 

N 

$ 


N 

N 

M 

Vi 


Ns 

«0 

N 

Vi 






A 


Qo 

N 

Os 

\ 




































































































2? 

«0 


































































































Q 


N 

1 


































































































u 


1 

fs 


































































































fQ 


1^ 




































































































N 

N 


































































































s. s 


0 






'n 

\ 


0 

N 


'o 

fNJ 


o 


Vi 




Vi 

V 






Vi 

vs 


'o 

\s 




Vi 

K 


o 


0 

Os 


0 

o 

N> 


0 

N 

n; 


0 

V 

n: 


O 

VQ 

N 


§ 

\ 


0 

N 





61 



TABLE IV: (Continued) 
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DISCUSSION OF RESULTS 



Comparing the power spectral solution and the theoretical solu- 
tion (Figs. 13, 14, 15, and 16) by matching the phase and amplitude 

A 

curves for -r— , it is noted that in both cases the phase shift of the 

ivy 

transfer function is in close agreement and the greatest amplitude of 

A 

•— , as established by both solutions, occurs at the same frequency, 
0.6 cycles per second, however a discrepancy in the value of the peak 

A 

amplitude of is apparent* The amplitude of the transfer function, 

as determined by the power spectral method, is consistently less, by 
an almost constant difference, than the amplitude as established by the 
theoretical solution for the same range of frequencies. The validity 
of the transfer function amplitude as established theoretically is con- 
sidered to be accurate Inasmuch as two different solutions produced 
identical results and the peak amplitude of 0.3 is a realistic value 
for an aircraft such as the F-100. Refer to Myslett and Smith (10). 

The power spectral solution, as derived, would imply that in flight 
the aircraft was not as disturbed by the turbulence as the theoreti- 
cal solution predicted. This is more obvious when examining the auto- 

« 

correlation function for Q as computed and graphed by the digital 
program. Although the distribution of the fluctuations of the auto- 
correlation function for Q was in agreement with theory; the magni- 
tude of the fluctuations was extremely small with an average absolute 
value of .000001. 

The computational results and the observation mentioned above 
led to a close scrutiny of the instrument sensitivity on the test air- 
craft used to gather the data. In particular, the sensitivities of 
the instruments measuring pitch angle, pitch rate, normal acceleration 
at the C.G., and angle of attack were checked because these parameters 
provide the most vital information regarding the actual dynamic response 
of the aircraft to turbulence. The entire flight record for Run number 
8 was re-examined for the typical inflight measurements of the para- 
meters mehtloned above. The results are summarized in Table V. 
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Table V: Instrument Sensitivity and In-Flight Measurements for 

Fit. 26, Run 8. 



Parameter Measured 


Inst. Sensitivity 
(in terms of trace 
deflection) 


Average Deflection 
Flight Record 


Typical 

Resultant 

Value 




12°/inch 


.22 inch 


t 1.6° 


o 


9.1° 

sec/ inch 


. 07 inch 


+ _o . 

- .7 /sec 




4g 's/inch 


.09 inch 


- ,36g 




4. 5°/ inch 


.15 inch 
.05 inch 


.2°(-o^»,) 



The tabulations above indicate a generally small response of the 
aircraft to the clear air turbulence encountered. The data from 
Table V, the diminished magnitude of the autocorrelation function for 
^ , and the consistent trend of the power spectral solution for the 

A ^ 

amplitude of to be less than the theoretical amplitude of 

indicate that t^e sensitivity of the instrumentation {board the test 
aircraft was insufficient. This conclusion was further verified during 
the investigation of the instrument sensitivity at NASA when discussion 
with the NASA engineers revealed that the aircraft used in collecting 
this data had been originally instrumented for the purpose of collecting 
badeuvering type data. The instrument sensitivity requirements for this 
type of maneuvering flight are not as stringent as those necessary for 
the power spectral type analysis of turbulence. The accuracy of the 
response is even more critical in a power spectral analysis since the 
associated digital computation is of such length that any initial dis- 
crepancies will be magnified by "rounding off" errors which are carried 
throughout the computer program. The amplitude discrepancy of the 
transfer function noted above, however, was of indirect value in an 
academic sense in that it generated a thorough secondary analysis of the 
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flight test data and a better appreciation of the entire flight test 
program from considerations of the Instrument location and sensitivity, 
to the digital program Itself. 
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SUMMARY AND CONCLUSIONS 



Despite the single discrepancy in the solution for the transfer 
function, 7— , it is felt that this analysis did represent a signi> 
ficant and unique application of power spectral methods because the 
transfer function between a random input and the resultant aircraft 
dynamic response was established on a continuous basis over a finite 
frequency range. The amplitude discrepancy between the theoretical 
and power spectral solutions for the transfer function is considered 
to be entirely due to improper instrument sensitivity and not to the 
mathematical computations or the power spectral analysis itself. A 
closer correlation between the amplitudes of the transfer functions 
would be realized using the same analytical procedures if the in- 
flight response measurements were more accurate. 

Future effort in turbulence response analysis, particularly with 
power spectral applications, will require that more sensitive instru- 
mentation be incorporated on the test aircraft. In addition, the angle 
of attack vane should be re-located from its present position on the 
vertical stabilizer to a position at least two and a half fuselage 
diameters ahead of the aircraft as proposed in Refs. 1 and 11. 

Although this analysis considered a singular application of 
the power spectral method, additional areas of investigation related to 
gust loading would be worthwhile for future research. A preliminary 
study, by Press and Houbolt (12), has been conducted to determine the 
variation in normal acceleration with center of gravity position for a 
fighter type aircraft. Through power spectral analysis, the ratio of 
the root mean square values of the normal acceleration, determined by 
the area under the gust spectrum, was readily established between center 
of gravity positions. It was determined that as the center of gravity 
position is shifted aft, the peak frequency response moves to a lower 
frequency and diminishes slightly. Other significant areas of future 
study and investigation include the effect of low short period damping 
on the gust loads of aircraft in continuous turbulence and the magni- 
fication in gust load stresses that results from wing flexibility. 

Refer to Press and Tukey (13). 
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In general, the power spectral methods of analysis may be re- 
garded as an effective and accurate tool for analysis of a random type 
disturbance. In particular, it is a powerful method for further ex- 
ploration of the science of gust loading. 
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Appendix A 



The Sequence of Digital Computer Steps in the Power Spectrum Analysis 



The mathematical steps outlined In Table VI explain the compu- 
tational procedure behind each of the digital computer steps In the 
power spectrum analysis (Ref. 14). All Information In Table VI refers 
to Fig. 17, and each step in the table Is related to the appropriate 
computer operation by the computer sequence number. 

The symbols used In Table VI are defined as follows: 

1^^ value of discrete, equl-spaced time series X(t) 

X^ 1^^ value of time series X(t) after subtracting mean 

1^^ value of time series Y(t) after subtracting mean 

P Lag 

n Number of discrete data points 

m Maximum lag 



(P) 

X 

(h) 

X 



SP 



xy 

,(h) 

xy 

,(h) 

-xy 

j(h) 

^xy 



SC 



(h) 

xy 



Constant time interval 

Autocovariance of series X at lag p, ■ R 

X X 



h n 
m z^t 



(h) 



Power spectral estimate of series X at frequency 

T,(h) „ h TT 

P ■ P nrzrt 

X X 

Smoothed power spectral estimate of series X at frequency 

h TT 
xa. ^ X. 

Cross-covariance between series X and Y at lag p 

Cross-power spectral estimate at frequency h 7T 

m A t 

Cospectrum, the real part of p^^^ at frequency h 7f 

m A t 

Quadrature spectrum, the Imaginary part of p^^^ at frequency 

xy 

h 77 
m A t 

Smoothed cospectrum at frequency h Tf 

m A t 
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Smoothed quadrature spectrum at frequency h ^ 



xy 



.(h) 



m 1 



am' ' Amplitude of cross spectrum at frequency h 7f 



xy 



m ^ t 



PHAS^*'^ Phase of cross-spectrum at frequency h 77' 

m A t 

tC*') Transfer function at frequency h 7 r 
xy 



m t 



,(h) 



TAM^y' Amplitude of transfer function at frequency h 



m 



7f 

t 
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Series / •OmEAN x- ^ AUTO RAW EST. ^ SMOOTHING ^ 

— CORRELATION ->■ OF ► POWER SPECTRUM r>- CHECKSUM 

(x=0) data function pwr.spec. by hamming 



IP 

53 h* 

^ <0 '' 


^TRANSFER 

FUNCTION 

PHASE 




1 i 


1 
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FIG. !7 SEQUENCE OF DIGITAL COMPUTER STEPS IN POWER SPECTRUM ANALYSIS 
/. X and y correspond to S and Wg respectively. 

2. Wg is computed from eqn.(5): Wg = V Ct^, -V 6 + Ja^ dt-H^O. 

3. Computer sequence numbers are circled. 








TABLE VI: Sequence of Digital Computer Steps In The Power Spectrum 

Analysis 



COIPUTER 

SEQUENCE 

NUMBER 


MATHEMATICAL OPERATION 


REFER 


TO PAGE: 


© 


The mean of the data Is first 
calculated and subtracted from each 
value of X' . 


(Original Data) 

e 76 




- /77^ 1= 




81 


(D 


The autocovariance is then computed « 

J- V 1/ 

= n-p 


o 

e 


77,78 








82,83 










® 


The raw estimate of the power spectrum 
Is obtained by 

(/,) y. ^ rp) /?p7T 

ZD _ l-€pR^ cos^^ 

^ 77 /Dro ^ ^ ^ 


o 

e 


79 




h//ye/^<s < p< m 




84 \ 












The raw estimates of the power spectrum 
are then smoothed by "hamming" 






! 


= . 34P^‘°^ ^ 


O 

6 


80- 


t 

\ 


5pf'=.Z5pJ'’'^+.5‘lfi’''^*.25/^‘'’'^o <h<m 






! i 

i ( 

k 


JPj"^= .54f^''”^ ^ 




85 
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TABLE VI (Continued) 



COMPUTER 

SEQUENCE 

NUMBER 


MATHEMATICAL OPERATION 


• 

REFER TO PAGE: 




A check sum is computed to check the 
computations o£ the estimates « 


h 79 




7T ft / (o) (f>)l 




This should equal R^^\ the autocovari- 
ance at zero lag. 


i 84 




The cross-covariance is computed by the 
formula 

0 


86,87, 

88 










The cross-spectrum, is given by 

/^xy - Cxy ^ Oxy 




© 

® : 


The cospectrum is obtained from 

The quadrature spectrum is obtained from 


89 

89 




^ (h)_ At y ^ (p 

Qxy rr fp^'\yy rX/y jSUl /-y-) 

p-o ' ' 






y, 2, ,f-n 




1 


ffz p=o,n-t 

€p / 

O < p<rn 




1 ® 

1 


Both the cospectrum and quadrature 
spectrum are smoothed by **hamming” as 
shown in step 4* 
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TABLE VI (Continued) 



C(»(PUTER 

SEQUENCE 

NUMBER 


MATHEMATICAL OPERATION 


REFER TO PAGE: 


@ 


The amplitude and phase of the cross- 
spectrum are computed. 


89,90 








@ 




89,91 




The phase angle computed is the phase 
shift of series Y with respect to 
series X. 




® 


The transfer function from series X 
to Y is given by 








92 


@ 


The amplitude is given by 






(h) (h) / (/,) 

rAM,y = 


92 




The phase shift from X to Y is the same 
as that of the cross-spectrum. Both 
the transfer function amplitudes from 
X to Y and Y to X are printed. 


i 
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6~HPD2T^lJT0CDyAftUNCE AND~TOHSR~ SPECTRXt^t TOm i S - VEKSIDW yiU» 7~ n; IW 

health sciences computing facility, UCLA 



~ ~ ~ P 1TD rTT“H~H-DrH I'B » RUK 1' 



INPUT DATA TO BE PRINTED OUT YES 






INPUT SERIES TO BE PLOTTED OUT ' 

OETRENOING » 


-t ~ ‘ 




PRBWMITENING 







VALW OrtONSTANT C“05EC» IN rHE~PPEbH1TEhING TRANSFORNATIOK ZITI-Xd/n^CXai >0. 

N UHBBft OF SERIES - 2 

NUMBER OF DATA POINTS • 920 

NUMBER OF LAGS CMGSE^N • 100 

NUMBER OF SELECTION CAROS « I 

USE PREVIOUS OATA 

COMStANT TIME interval • 0. 10000 second 

NUMBER OF VARIABLE FORMAT CARDS ■ I 



1 






original data of series 



- 0.00061 


- 0. 00066 


0.00027 


0.00027 


0.00027 


0*00027 


< 0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00002 


0.00002 


0 . d 0002 


0 .0000 2 


0.00002 


0.00002 


0.00002 


0.00002 


0.00052 


0.00102 


0.00202 


0.00290 


0.00315 


0.00260 


0.00152 


0.00066 


0.00001 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00023 


0*00046 


0.00190 


O.OOJIS 


0.00603 


0.00378 


0.00278 


0.00160 


0.00077 


0.00016 


- 0.00023 


- 0.00111 


- 0.00211 


- 0.00211 


- 0.00211 


- 0.00169 


- 0.00068 


- 0.00023 


0.00115 


0.00115 


0.00115 


0.00190 


0.00177 


0.00090 


0.00002 


- O.OOOU 


- 0.00011 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00089 


- 0.00136 


- 0.00269 


- 0.00269 


- 0.00269 


- 0.00068 


0.00066 


0.00127 


0.00265 


0.00202 


0.00152 


0.00077 


0.00039 


- 0.00023 


- 0.00023 


- 0.00024 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00023 


- 0. 00023 


0.00052 


0.00002 


- 0.00099 


- 0 . 00 O 84 


- 0.00023 


- 0. 00073 


- 0.00086 


- 0.00161 


- 0.00231 


- 0.00199 


- 0.00199 


- 0.00161 


- 0 .00136 


- O.OOOU 


0.00260 


0.00628 


0.00278 


0.00278 


0.00215 


0.00160 


0.00039 


0.00039 


0.00039 


0.00066 


0.00115 


0.00152 


0.00152 


0.00152 


0.00360 


0.00616 


0.00614 


0.00056 


0.00656 


0.00628 


0.00165 


- 0.00011 


- 0.00161 


- 0.00369 


- 0.00387 


- 0.00387 


- 0*00362 


- 0.00036 


0.00039 


0.00177 


0. 00603 


0*00516 


0.00667 


0.00667 


0.00606 


0.00253 


- 0.00068 


- 0.00211 


- 0.00637 


- 0.00637 


- 0.00637 


- 0.00269 


- 0.00126 


C . 00016 


0.00016 


0.00016 


0.00016 


0.00016 


- 0.00061 


- O.OOOU 


- 0.00036 


- 0.00036 


- 0.00036 


- O . C 0036 


- 0 . C 0086 


- 0.00036 


- 0.00036 


- 0.00036 


- 0.00034 


0.00077 


0.00616 


0.00729 


0.00662 


0.00630 


0.00830 


0.00729 


0.00603 


0.00152 


0.00016 


- 0.00084 


0.00002 


0.00090 


0.00090 


0 . C 0090 


0.00090 


0.00066 


0.00016 


0.00016 


0.00016 


0. 00016 


0.00016 


0.00066 


0.00115 


0.00115 


0.00102 


0.00052 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


0.00027 


- 0.00073 


- 0.00111 


- 0.00161 


- 0.00186 


- 0.00164 


- 0.00099 


- 0.00023 


0 . 0 C 016 


0.00066 


0.00039 


- 0.00036 


- 0.00111 


- 0.00186 


- 0.00269 


- 0.00311 


- 0.00312 


- 0.00237 


- 0.00111 


- 0 . 00 C 41 


- C . 00011 


- O.COOll 


- 0.00011 


- 0.00011 


- o.oooii 


- 0.00068 


- 0.00086 


- 0.00066 


- 0.00086 


- 0.00066 


- 0.00066 


- 0.00086 


- 0.00136 


- 0.00199 


- 0.00226 


- 0.00226 


- 0.00199 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00161 


- 0.00211 


- 0.00369 


- 0.00362 


- C . 00312 


- 0.00326 


- 0.00687 


- 0. 00312 


- O.OOiU 


- 0.00126 


0.00190 


0. 00661 


0.00691 


0.00492 


0.00629 


O . O 031 S 


0.00039 


- 0.00061 


- 0.00126 


- 0.00111 


- 0.00111 


- 0.00169 


- 0.00211 


- 0.00214 


- 0.00199 


- 0.00161 


- 0.00099 


- 0.00036 


- O.COOll 


0.00052 


0.00090 


0.00090 


- 0.00011 


- O.OOIU 


- 0.00161 


- 0.00262 


- 0.00269 


- 0.00186 


- 0.00099 


- 0.00068 


0.00016 


0.00027 


0.00016 


- 0.00014 


- 0.00068 


- 0.00066 


- 0.00136 


- 0.00211 


- 0.00126 


- 0.00086 


- 0.00099 


- 0.00099 


- 0. 00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00099 


- 0.00098 


- 0.00086 


0.00016 


0.00039 


0.00039 


0.00039 


0.00027 


0.00002 


- 0.00086 


- 0. 00199 


- 0.00267 


- 0.00287 


- 0.00267 


- 0.00287 


- 0.00269 


- 0.00186 


- 0.00136 


- 0.00111 


- 0.00111 


- 0.00136 


- 0.00164 


- 0.00199 


- 0.00199 


- 0.00199 


- 0.00199 


- 0.00186 


- 0.00161 


0.00016 


0.00090 


0.00260 


0.00391 


0.00603 


0.00364 


0. 00177 


- 0.00036 


- 0.00169 


- 0.00161 


- 0.00287 


- 0.00287 


- 0. 00287 


- 0.00287 


- 0.00237 


- 0.00161 


- 0 . 0 C 086 


- 0.00061 


- 0.00036 


- 0.00023 


0.00016 


0.00102 


0.00090 


0*00052 


- O.OOOU 


- 0.00023 


- 0.00023 


- 0.00023 


- 0.00023 


- 0*00099 


- 0.00161 


- 0.00161 


- 0. 00169 


- 0 . 0 Q 16 S 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00169 


- 0.00226 


- 0.00625 


- 0. 00563 


- 0.00563 


- 0.00687 


- 0.00337 


- 0.00184 


- 0.00136 


- 0.00066 


- 0.00086 


- 0.00061 


- 0.00044 


- q^floouT 


0.00039 


0.00115 


0.00190 


0.00177 


0.00177 


0.00027 


0.00002 


0.00002 


0.00002 


O . 00 O 90 


0.00253 


0.00278 


0.00353 


0.00360 


0.00360 


0.00366 


0.00228 


0.00066 


0.00002 


- 0.00068 


- 0.00061 


- 0.00061 


- 0.00066 


- 0.00068 


- 0.00068 


- 0.00068 


- 0.00068 


- 0. 00068 


- 0.00068 




0. 00115 


0. 00260 


0.00616 


0.00616 


0.00190 


- 0.00011 


- 0.00073 


- 0.00226 


- 0.00195 


- 6 > 00 M 9 


- 0.00126 


- 0.00099 


- 0.00023 


0.00077 


0.00077 


0.00027 


- 0.00086 


- 0.00073 


- 0.00074 


- 0. 00073 


- 0.00176 


- 0.00176 


- 0.00176 


- 0.00176 


- 0.00126 


- 0.00073 


- 0 . O 0036 


- 0.00023 


- 0.00024 


- 0.00023 


- 0.00099 


- 0.00099 


- 0.00111 


- 0.00169 


- O.OOlll 


- 0.00068 


0.00002 


0.00090 


0.00127 


0.00127 


0. 00090 


0.00052 


0.00052 


0 . C 0052 


O . C 0016 


- 0.00068 


- 0.00066 


- 0.00126 


- 0.00126 


- O.OOlll 


- 0. 00061 


0.00052 


0 . 00 CS 2 


0 . 000 S 2 


b . maisi : 


0.00016 


“^ o . oooiiT ” 


- 0.00136 


- 0.00184 


"'-(I i flow 4 


- 0.00061 


0.00115 


0.00290 


0. 00290 


0*00278 


0.00202 


0.00052 


- 0.00073 


^ 0.00326 


- 0.00675 


- 0.00586 


- O .00525 


- 0.00276 


- O . CO 086 


0.00052 


0.00077 


0.00265 


0.00628 


0.00678 


0.00661 


0.00360 


0.00177 


0.00039 


- 0.00011 


- 0*00126 


- 0.00226 


- 0.00211 


- 0.00211 


- 0 . 0 Q 2 U 


- 0.00111 


- 0. 00136 


- 0.00061 


O . OO 039 


0.00160 


0.00315 


0.00653 


0.00653 


0. 00653 


0.00651 


0.00653 


0.00391 


0.00165 


0.00016 


0.00016 


- 0.00034 


- 0.00023 


0*00090 




0*00678 


QtOQ?6i 



0.00566 

- 0.00675 

0.00002 



0.00321 

- O . O 0663 

0.00202 

- 0.00036 

- 0.00126 

0.00190 

- 0 . 0006 B 

0.00102 

0.00052 

- 0.00161 

0.00260 

- 0.00262 
0. 00260 
- 0.00166 
0.00016 
- 0 . C 0 O 66 
- 0.00169 
0.00376 
0.00077 
- 0.00011 
0.00016 
- 0.00073 



0.00691 
- 0.00186 
0. 00016 



0.00653 
- 0. 00575 
0.00052 
O ^. OcO^t 

- 0.00169 
0.00115 
' o.OOUS 
- 0.00068 
0.00127 
6.00177 
- 0*00099 
- 0. 00011 
- 0.00061 
- 0.00136 
0.00066 
-0.00186 
0.00002 
- 0.00023 
-0.O0161 
- 0.00011 
0.00077 
- 0*00023 
- 0.00011 
- 0.00086 



0.00315 

0 . O 0016 

0. 00016 

- 0.00237 
-^.000 11 
0.00303 
- 0.00169 
0.00027 

TToom 
- 0.00068 
0.0 0160 
o . OoI ^ 
- 0.00023 
- 0.0 0169 
— 0. 0006 1 
- 0.00073 
- 0.00011 
^. O 0 i 86 ' 
- 0.00036 
- 0.00011 
-0.00176 
- 0.00023 
0.00127 
- O . 0 O ^ Cl 6 
- 0.00023 
- 0.00136 



0.00090 
C . 00090 
0.00016 
■“ OTucm - 
- 0.00126 
- 0.00086 
o.ffojTr 

- 0.00169 

- 0.00073 

- 0.00068 
0. 00102 _ 
■ JT . OOITb 
0.00077 
- 0 .0037 5 

- o.ooiir 

- 0.00036 
- 0.00061 
“-0. 00126 
- 0.00086 
0.00090 
- 0.00186 
- 0.00036 
0.00190 
- G . 00 G 36 
- 0.00036 
- 0.00111 



- 0.00011 
0.00260 
0.00016 
' 0.00328 
- 0. 00068 
- 0.00267 
?T.0025r 
- 0.00169 
- 0.00073 



- 0.00136 

0.00260 

0.00016 

~DT07T25T" 

0.00039 

- 0.00287 

O . O 0228 

- 0.00073 

- 0.00073 



- 0.00068 - 0.00068 

0 .00027 - 0. 00036 

( T .~ Q 0378 0.00278 

0.00228 0.00378 

- 0.00337 - 0. 00299 

=^. U 026'2 - 0:00262 

0.00177 0.00391 

- 0. 00099 - 0.00176 



- 0. 00369 
0.00260 
0.00016 
6.0C086 
0.00165 
- 0.00267 
0.00177 
- 0.00068 
- 0.00068 

- o ! o 0036 
- 0 .00073 
0^.00102 
0.00579 
- 0. 00237 
- 0.00262 



- 0.00666 
0.00260 
- 0.00036 
- 0.00199 
0. 00260 
- 0.00287 
0.00077 



- 0.00066 
0.00160 
-^OOITF 
- 0. 00036 
0.00152 



- 0.00036 

- 0.00061 



- 0.00066 

0.00115 



-0.00073 
- 0.00061 
0.00090 
0. 00016 
- 0.00036 
- 0.000 23 



0.00691 
- 0.0017 6 
- 0 . 0006 T 
- 0.00099 
0.00039 
0.00152 
- 0.00061 
0.00066 
0.00077 
- 0.00068 
- 0. 00023 



- 0.00023 

- 0.00023 

- 0.00036 

- 0.00176 

0.00729 
- 0. 00176 
-0 .00262 
0.00691 
- 0.00166 

- 0.00099 

- 0.00011 

0.00178“ 

- 0.00036 

0.00016 

0.00077 

- 0.00061 

- 0.00036 



- 0.00639 

0.00202 

- 0.00011 

- 0.00387 

0.00265 

- 0.00136 

0.00059 



0.00077 

0.00016 



- 0.00023 
- 0. 00126 



0.00656 

- 0.00073 



0.00691 
- 0. 00166 



- 0.00099 

- 0.00073 



0.00391 

- 0.00036 

- 0.00011 

- 0.00061 

- 0.00066 



- 0.00625 
0.00077 
0.00165 
” 0 . C 04 I 3 “ 
0.00265 
- 0.00066 
-0. 00048 
0.00177 
0.00077 

0.00027 

- 0 .00061 

0I0OSO6 
- 0 .00068 
1.00262 
0. 00661 
- 0 .00166 

- 0.00066 

- 0.00066 



0.00391“ 

0.00027 

- 0.00011 

- 0.00061 

- 0.00111 



- 0.00141 


-0.00176 


- 0.00176 


- 0.0 0 176 


- 0.00161 


-0 «00041 


0 . 0009 Q 


0.00260 


0.00603 


O . G 0603 


0.00391 


0.00265 


0.00160 


- 0.00036 


- 0.00066 


- 0.00176 


- 0.00237 


- 0. 00269 


- 0.00262 


- 0.00262 


- 0.00211 


- 0.00136 


- 0.00061 


- 0.00016 


0.00016 


0.00066 


0. 00160 


0.00290 


0.00315 


0.00290 


0 . 0^15 


0. 00346 


O . OO 260 




O . O0O02 


- 0.00011 


- 0.00011 


- O.OOOU 


- 0.00011 


- O.OOOU 


- 0.00068 


- 0.00099 


- O . OOC 86 


- 0.00099 


- 0.00099 
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graph of AUT0C0VARIANC6 FUNCTICN OF SERIES 2 PLOTTEG AGAINST TIME UP TO A LAG OF 10.0000 SECOND 
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6.600 


, 






4 




6.600 


6.700 








4 




6.700 


6.800 


• 






4 




6.800 


6.900 








4 




6.900 


7.000 


4- 






4 




7.000 


7.100 

7.200 


• 






4 

4 


• 


7.100 

7.200 


7.300 


. 






4 




7.300 


7.400 


, 






4 




7.400 


7.500 


♦ 






4 




7.500 


7.600 


, 






4 




7.600 


7.700 


. 






4 




7.700 


7.800 


, 






4 




7.800 


7.900 








4 




7.900 


8.000 


♦ 






4 


■f 


8.000 


8.100 


• 






4 




8.100 


8.200 


, 






4 




8.200 


8.300 


. 






* 




8. 300 


8.400 


. 






4 




8.400 


8.500 


♦ 






4 


♦ 


8.500 


8.600 


• 






4 




8.600 


8.700 


• 






4 




8.700 


8.800 


• 






4 




8.800 


8.900 


, 






♦ 




8.900 


9.000 


♦ 






« 


♦ 


9.000 


9. 100 


, 






4 




9.100 


9.200 


. 






4 




9.200 


9.300 








* 




9.300 


9.400 


, 






4 




9.400 


9.500 

9.600 


♦ 


— 




4 

4 




9.500 
9. 600" 


9.700 


, 






4 




9. JOO 


9.800 


, 






4 




9. BOO 


9.900 


. 






4 




9. 900 



- 0.006 



-0.00t> 



- 0.000 











<g) 


FREQUENCY 


CO-SyfCtRUH 


QUAORAIURE SPECTRUM 


“AN<>lITUDE of CR«S-SPEdTRUH 


PHASE OF CROSS-SPECTRUM 


ICYCLES/SECONOI 


OF 


OF 


OF 


OF 




SERIES 1 ANO 2 


SERIES 1 AND 2 


SERIES 1 ANO 2 


SERIES 1 ANO 2 * 


-0. 


-0.0021232 


-0.0032847 


0.0039111 


0.6586698 


0.050 


-0.001A642 


-0.0036339 


0.0039176 


0.6890397 


0. 100 


-0.0007291 


-0.0016674 


0.0020233 


0.6913315 


0.150 


-0.0002056 


-0.0016613 


0.0016740 


0.7303643 


0.200 


0.0000641 


-0.0015001 


0.0015024 


0.7569146 


0.250 


0.0001300 


-0.0013435 


0.0013498 


0.7653466 


0.300 


O.OC 00567 


-0.0009606 


0*0009624 


0.7591837 


0.350 


0.0001166 


-0.0010761 


0*0010626 


0.7674672 


0.400 


0.0001705 


-0.0006342 


0.0006567 


0.7917943 


0.450 


0.0002059 


-0.0002335 


0.0003113 


0.6650162 


0.500 


0.0002611 


-0.0000399 


0.0002639 


0.9775740 


0.550 


0.0002331 


0.0001256 


0.0002649 


0.0787852 


0.600 


0.0003949 


0.0006309 


0.0007443 


0.1609962 


0.650 


0.0006710 


0.0006557 


0.0010674 


0.1441602 


0.700 


0.0004431 


0.0006563 


0*0007916 


0.1554665 


0. 750 


0.0001671 


0.0004414 


0.0004794 


0.1861661 


0.600 


0.0000596 


0.0004297 


0.0004336 


0..2260736 


0.650 


0.0001251 


0.0003221 


0.0003455 


0.1910421 


0.900 


0.0001735 


0.0002061 


0.0002709 


0.1393655 


0.950 


0.0002002 


0.0001729 


0.0002645 


0. 1133860 


1.000 


0.0001696 


0.0002451 


0.0003100 


0,1451075 


1.050 


0.0001140 


0.0001750 


0.0002068 


0*1560855 


1.100 


0.0000712 


0.0001851 


0.0001963 


0. 1915373 


1.150 


0.0000592 


0.0001267 


0.0001417 


0.1613452 


1.200 


0.0000335 


0.0000619 


0.0000665 


0.1662467 


1.250 


0.0000442 


0.0000451 


0.0000632 


0,1265222 


1.300 


0.0000756 


0.0000637 


0.0001129 


0.1329414 


1.350 


0.0000363 


0.0000365 


0.0000543 


0.1254156 


1.400 


0.0000010 


0.0000441 


0.0000441 


0.2465460 


1.450 


0.0000006 


0.0000140 


0.0000141 


0.2411736 


1.500 


0.0000329 


0.0000306 


0.0000449 


0.1190123 


1.550 


0.0000403 


0.0000060 


0.000041 1 


0.0311476 


1.600 


0.0000137 


0.0000216 


0.0000256 


0*1605606 


1.650 


0.0000146 


0.0000063 


0.0000159 


0.0647518 


1.700 


0.0000304 


0.0000216 


0.0000375 


0.0991162 


1.750 


0.0000267 


0.0000012 


' 0.0000267 


0.0073651 


1.600 


0.0000251 


0.0000235 


0.0000344 


0.1197049 


1.650 


0.0000391 


0.0000061 


0.0000399 


0.0323369 


1.900 


0.0000310 


0.0000201 


0.0000370 


0.0915918 


1.950 


0.0000136 


-0.0000020 


0.0000139 


0.9776278 


2.000 


0.0000067 


0.0000106 


0.0000139 


0.1416096 


2.050 


-0.0000002 


-0.0000041 


0*0000041 


0*7435969 


2.100 


0.0000019 


0.0000169 


0.0000 170 


0*2321246 


2.150 


0.0000097 


0.0000053 




0.0^9590 




b. ooooio'/ 


070*00 50‘9iT 




*"'“*^* ^0S’doi45 


on rff6'76fi~ ^ 


2.250 


0.0000076 


-0.0000064 


0-0000099 


0.6664183 


2.300 


0.0000043 


0.0000063 


0.0000093 


0.1733696 


2.350 


0.0000069 


-0.0000064 


0.0000094 


0.6614613 


2.400 


0.0000071 


0.0000036 


0.0000061 


0.0779053 


2.450 


0.0000066 


-0.0000047 


0.0000061 


0.9016696 


2.500 


0.0000046 


0.0000092 


0.0000104 


0*1726165 


2.550 


0.0000041 


-0.0000039 


0.0000056 


0.6762449 


2.600 


0.0000033 


0.0000053 


0.0000063 


0.1610474 


2.650 


0.0000072 


-0.0000031 


0.0000076 


0.9357431 


2.700 


0.0000077 


0.0000062 


0.0000099 


0*1062167 


2.750 


0.0000014 


-0.0000033 


0.0000036 


0.6121454 


2.600 


0.0000032 


0.0000070 


6.0000077 


0*1811776 


2.650 


0.0000147 


-0.0000011 


0.0000147 


0.9664765 


2.900 


0.0000170 


0.0000034 


0.0000 17 3 


0.0315632 


2.950 


0.0Q00060 


-0.0000116 


0.0000142 


0.6446193 


3.000 


0.0000007 


-0. 0000063 


0.0000063 


0.7666611 


3.050 


-0.0000049 


-9*0000076 


0.0000093 


0.6603565 


3.100 


-0.0000064 


0.0000040 


0.0000076 


0.4116942 


3.150 


-0.0000052 


-0.0000037 


0.0000064 


0*5991808 


3.200 


0.0000003 


0.0000006 


0.0000007 


0.1898396 


3.250 


b. 00000 13 


-0. 0000070 


0 . OOCTDOn 


“077700337 


3.300 


0.0000014 


0.0000029 


0.0000032 


0.1772105 


3.350 


0.0000050 


0.0000003 


0.0000050 


0.0097569 


3.400 


0.0000046 


0.0000053 


0.0000072 


71.TTn«75 


3.450 


0.0000024 


-0.0000042 


0.0000046 


0.6324596 


3.500 


0.0000003 


0.0000009 


0.0000010 


0.2016066 


3.550 


-0.0000004 


-0. 0000044 


0.0000045 


0.7342104 


3.600 


-0.0000026 


-0.0000002 


0.0000026 


0.5132156 


3.650 


-0.0000032 


-0.0000050 


0.0000059 


0.6595750 


3.700 


-0.0000022 


0.0000014 


0.0000026 


0.4075344 


3.750 


-0.0000005 


-0.0000026 


0.0000026 


0.7229510 


3.600 


0.0000007 


-0.0000001 


0.0C00007 


0.9710475 


3.650 


0. 13000005 


-0.0000067 


0.0000067 


O.T624504 


3.900 


-0. 0000003 


-0.0000031 


0.0000031 


0.7356646 


3.950 


-0.0000006 


-0.0000036 


0.0000037 


C. 7246844 


4.000 


0.0000015 


-0.0000006 


0.0000016 


0.9367036 


4.050 


-0.0000004 


-0.0000036 


0.0000036 


0. 7319391 


4.100 


-O.OCO0O24 


0.0000010 


0.0000026 


0.4360235 


4.150 


-0.0000014 


-0.0000042 


0.0000044 


0.6974187 


4.200 


-0.00C0006 


-0.0000040 


0.0000041 


0.7273278 


4.250 


0.C000004 


-0,0000050 


0.0000050 


C. 7632795 


4. 300 


O.OOOC014 


0.0000002 


0.0000014 


C. 0195076 


4. 350 


O.CCCOC06 


-C. 0000025 


0.0000026 


C. 7665616 


4.400 


O.OC00016 


-0.0000016 


0.0000024 


0.6642746 


4.45C 


0.CCCCC21 


-0.0000049 


0.CC00053 


C. 8155177 


4. 500 


C.CC0CO15 


-0.C000020 


O.OC0OO2S 


C. 8491764 


1.550 


O.OCOCOOb 


-0.0000007 


O.OCOOOlO 


0-8694575 


4.600 


C.00OC016 


-0.0000001 


0.0000016 


0.9939513 


4.650 


O.CCOC016 


-0.0000006 


0.0C00016 


0.9235646 


4.700 


c.ccccccc 


0.0000025 


O.CC00025 


C. 24748?? 


4.7^0 


C.CC0C0C2 


O.COOOOT1& 


O.OXXOTTOr 


"170156656“ 


4.600 


C.CCCCCC6 


-O.COOO0O6 


O.OCOOOlO 


C.84637C9 


•..650 


0.CCCCC24 


0.0000001 


O.OOOCC24 


O.0C4456C 


-.90u 


C.00U0o?5 


0.0000014 


0.0000026 


o. 06154 7fc 


. 0 


u.0000';2 ) 


-0.0000006 


0.00000^ 3 


0.9586746 




> .CCC(A 


-0.0000015 


O.CCCOC 3C 


«.,9 144 ‘ f t 



L' 



IRC. r 



GRAPH OF AMPLlTUOe OF CROSS-SPECTRUM OF SERIES I AND 2 UN LOG SCALE) AGAINST FREOIENCV (CYCLES/SECONO ) 



-15.598 _ -13.587 -11.576 -9.565 -7.553 -5.542 

-14.593 -12.581 -10.570 -8.559 -6.548 
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PHASE OF CROSS'SPECTRUH OF SERIES 1 ANO 2 — PLOTTEO IN FRACTIONS OF A CIRCLE AGAINST FREQUENCY tCYCLES /SECONO) 

1*^00 



-0. ♦ 

0.023 . 
0.030 . 
0.073 . 
0.100 . 
0.123 ♦ 
0.130 . 
0.173 . 
0.200 . 
0.223 . 
0.230 ♦ 
0,273 . 
0.300 . 
0.323 . 
0.330 . 
0.373 ♦ 
0.400 . 
0.423 . 
0.430 . 
0.473 . 
0.300 ♦ 
0.323 . 
0.330 . 
0.373 . 
0.600 . 
0.623 ♦ 
0.630 . 
0,673 . 
0.700 . 
0.723 . 
0.730 ♦ 
0.773 . 
0.800 . 
0.823 . 
0.830 . 
0.873 ♦ 
o'. 900 . 
0.925 . 
0.930 . 
0,973 , 

1.000 ♦ 
1.023 . 
1.030 . 

1.100 . 
1.123 ♦ 
1.150 . 
1.173 . 
1.200 . 
1.223 . 
1.230 ♦ 
1.273 . 
1.300 . 
1.323 . 
1.330 . 
1.373 ♦ 
1.400 . 
1.423 . 
1.430 . 
1.475 . 
1.500 ♦ 
1.525 . 
1.330 . 
1.575 . 
1.600 . 
1.625 ♦ 
1.650 . 
1.675 . 
1.700 . 
1.725 . 
1.750 ♦ 
1.775 . 
1.800 . 
1.823 . 
1.830 . 
1.875 ♦ 
1.900 . 
1.925 . 
1.950 . 
1.975 . 

2.000 ♦ 
2.025 . 



-0.300 

+ 10 



0. 300 



AjlIQQ. 



1.300 






- 0 ^ 
0.023 
0.030 
0.075 
0.100 
0. 125 
0,150 
0.175 
0.2 00 
0.225 
0.230 
0.273 
0.30^ 
0.325 
0.330 
0.375 
0.400 
0.423 
0.430 

6.475 
0.300 
0.325 
0.550 
0.575 
0.600 

0.630 

0.675 

0.700 

0.725 

0.730 

0.775 

0.800 

0.823 

0.050 

0.873 

0.900 

0.975 

0.930 

0. 975 
l.OOIT 

1.025 
1.030 

IV075-- 

1.100 

1.125 

1.150 

1.175 

1.200 

1.225 

1.230 

1.275 

1.300 

1.325 

1. 350 
1.375 
1.400 
1.425 
1.430 

1.475 
1.300 
1.525 
1.330 
1.375 
1.600 
1.625 
1.630 
1.675 
1.700 
1.723 
1. 730 
1.775 
1.800 
1.823 
1.830 
1.875 
1.900 
1.925 
1.950 
1.975 
2.000 

2.025 
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0) 






FREQUENCY 


COHERENCE SQUARE 


AHPLITUOE OF 


AMPLITUOE OF 


PHASE OF 




ICrCLES/SECONOI 


OF 


TRANSFER FUNCTION 


TRANSFER FUNCTION 


transfer function 


^t/>vcr/oV 




SERIES 1 Ami 


^RQM 1 TQ 1 


from Z TQ 1 


FROM 1 TQ. 2 ..* 


i TO Z, 


“0. 


0.0619165 


6074. 7177734 


0.0000102 


0.6586698 


/Z3 O 


0.0500 


0.1614U7 


7526.7287598 


0.0000214 


0.6890397 


// 2.J 


0. 1000 


0.5300088 


4609.3800659 


0.0001150 


0.6913315 


///. / 


0.1500 


0.6889446 


3483.0234375 


0.0001978 


0.7303843 


97 / 


0-2000 


0.7242599 


2872,0730286 


0.0002522 


0.7589148 


aa 9 


0.2500 


0.9594974 


2880.1491394 


0.0003331 


0.7653466 


d^. 7 


0.3000 


0.6903551 


1847.2368774 


0.0003737 


0.7591837 


i. 


0.3500 


0.6757946 


1483.8199463 


0.0004554 


0.7674672 


S3, a 


0.^000 


0.4400037 


1052.9419556 


0.0004179 


0.7917943 


75.0 


0.^500 


0.2456576 


733.9263306 


0.0003347 


0.0650182 




0.5000 


0.2757262 


719.2819901 


0.0003033 


0.9775740 


__L 


0.5500 


0. 1686972 


391.6125259 


0.0004308 


0.0787852 


Zi 


0.6000 


0.5212792 


454.7567940 


0.0011463 


0.1609982 


52.5 — ▼ 


0.6500 


0.6685881 


5U.1065Z65 


0.0013081 


0. 1441602 


5/ 9 


0.7000 


0.6822397 


648.J1295166 


0.0010526 


•0. 1554885 


ss 9 


0.7500 


0.6187247 


978.2851639 


0.0006325 


0. 1861861 


C7 0 


0.8000 


0.6741928 


1280.6452759 


0.0005264 


0.2280736 


82. Z 


0.8500 


0.4412994 


1009.3552246 


0. 0004372 


0.1910421 


C8.8 


0.9000 


0.5377295 


1049.2672119 


0.0005125 


0.1393655 


50.2 


0.9500 


0.5397015 


1221.4488831 


0.0004418 


0.1133860 


4o.f 


1.0000 


0.6581574 


1486.2371521 


0.0004428 


0.1451075 


S2.3 


1.0500 


0.5283800 


1569.1015778 


0.0003367 


0.1580855 


5C. 9 


1.1000 


0.6613693 


1821.7242432 


0.000 3630 


0.1915373 


C8.9 


1.1500 


0.4437273 


1940.5514069 


0.0002287 


0. 1813452 


£5.3 


1.2000 


0. 6364876 


2976.3462524 


0.0002138 


0.1882467 


a 7. 7 


1.2500 


0. 3351934 


1872.0773163 


0.0001790 


0.1265222 


45.5 


1.3000 


0.5603038 


2384.6016235 


0.0002350 


0.1329414 


47.8 


1.3500 


0. 1766725 


1400.0268097 


0.0001262 


0.1254156 


4S / 


1.^000 


0.35671L1 


1176 a 2414616 


Q.QQQI8Q5 


<L.2.4.fe5-^8Q 


98.7 


1.4500 


0.0603871 


886.6184082 


0.0000681 


0.2411736 


S£. 8 


1.5000 


0.3647295 


2183.8581543 


0.0001670 


0.1190123 


42.5 


1.5500 


0.3726450 


2263.4462585 


0.0001646 


0.0311476 


// 2 


1.6000 


0.335729r 


1572.8514557 


UTJJUUTm 


0.1605806 


S7.t 


1.6500 


0.0976866 


930.6552048 


0.0001050 


0.0647518 


23 3 


1.7000 


0.5087685 


2947.2337952 


0.0001726 


0.0991182 


JS7 


1.7500 


0.3973063 


3354.8771667 


0.0001186 


070073651 


2.7 


1.8000 


0.5725720 


3776.2563477 


0.0001516 


0.1197049 


48. f 


1.8500 


0.5433132 


2517.9406433 


0.0002158 


0.0323389 


Ji.£ 


1.9000 


0.6925161 


2(526.1054230 


0.0003418“ 


07WT59TB 


33.0 


1.9500 


0.2104636 


1324.4456482 


0.0001589 


0.9776278 


JS/.t 


2.0000 


0.4060428 


3108.2918091 


0.0001306 


0. 1418096 


S/. o 


2.0500 


0.0568191 


1423.4937134 


070000109“ 


0.7635969 




2.1000 


0. 4143999 


2534.3022461 


0.0001635 


0.2321248 




2.1500 


0.1160870 


1467.5696869 


0.0000791 


0.0799590 




2.20ffiy 


0. 3373034 


JOlU. 008 1787 


07000TT7T 


0. 1 1 86766 




2.2500 


0.339L165 


2542. 1031189 


0.0001334 


0.8884183 




2.3000 


0. 3951004 


2339.2106323 


0.0001689 


0.1733896 




2.3500 


0.2073879 


1897.1470947 


0. 0001093 


“ 0^.88146n “^ 




2.4000 


0.0762371 


1174. 1041107 


0.0000649 


0.0779053 




2.4500 


0.0967490 


1012. 1558762 


0.0000956 


0.9016696 




2.5000 


0.5335481 


1869.1534729 


0.0001764 “ 


0:i728165 




2.5500 


0. 1620380 


1876.8699493 


0.0000863 


0.0782449 




2.6000 


0.3285159 


2753.6622009 


0.0001193 


0.1610474 




2.6500 


o.wms 


3r3B. 3084717 


0.0001319 ■ 


0.9T57431 




2.7000 


0.5032876 


4148.7296143 


0.0001213 


0.1082187 




2.7500 


0.0605351 


1708.0550995 


0.0000354 


0.8121454 




2.8000 


0. 1319702 


205ffi92l5393 


0. 0000641 ' 


0.1811776 




2.8500 


0.3315270 


3962.8581238 


0.0000837 


0.9884785 




2.9000 


0.4106817 


3165.8064880 


0.0001297 


0.0315632 




2.9500 


0.3886^681 


T?3^.272765T 


0.0001727 


0.8448193 


— 


3.0000 


0.0876520 


1178.8600616 


0.0000743 


0.7666611 




3.0500 


0.2599313 


1767.2790833 


0.0001471 


0.6603585 




3.1000 


0. 20671T8 


158^.5707428 


0.0001297 


0.4118942 




3. 1500 


0.1158226 


1444.4833832 


0.0000802 


0.5991808 




3.2000 


0.0019717 


197.2228985 


0.0000100 


0.1898396 




3.2500 


0.4477476 


3288.1257629 


0.0001362 


0.7788337 




3.3000 


0. 1373497 


1917.8540192 


0.0000716 


0.1772105 




3.3500 


0.2815928 


2685.6484680 


0.0001048 


0.0097569 




3.4000 


0.6578453 


2743.4693298 


<).0002398 


6.1331675 




3.4500 


0.5462859 


2131.6284485 


0.0002563 


0.8324596 




3.5000 


0.0266906 


534.1774063 


0.0000500 


0.2018086 




3.5500 


0.3734227 


2274.2352905 


0.0001642 


0.7342104 




3.6000 


0.0917691 


1079.4534607 


0.0000850 


0.5132156 




3.6500 


0.4369779 


2512.2828369 


0.0001739 


0.6595750 




3.7000 


0.0951813 


1324.5124817 


0.0000719 


0.4075344 




3.7500 


0.1409936 


2000.5970917 


0.0000705 


0.7229510 




3.8000 


0.0124684 


623.6468735 


0.0000200 


0.9710475 




3.8500 


1.0855527 


6093.4472656 


0.0001782 


0.7624584 




3.9000 


0. 1534083 


2117. 1085510 


0.0000725 


0.7358648 




3.9500 


0.2499397 


2859.4440918 


0.0000874 


0.7248844 




4.0000 


0.0353635 


935.8715439 


0.0000378 


0.9387038 




4.0500 


0. 3604958 


2818.8566284 


0. 0001279 


0.7319391 




4. 1000 


0.1919961 


1806.6541748 


0.0001063 


0.4360235 




4.1500 


0.2832790 


1950.9380035 


0.0001452 


0.6974187 




4.2000 


0.1815312 


1961.1654358 


0.0000926 


0.7273278 




4.2500 


0.5055732 


5117.8895264 


0. 0000988 


0.7632795 




<•.3000 


0.0627275 


1609.4082031 


0.015155m 


0. '6195076 




4.3500 


0.2640604 


2948. 5805054 


0.0000895 


0.7865818 




4.4000 


0. 1507886 


2235.9115295 


0.0000674 


0.8842748 




4. 4500 


0.5587885 


3309.3141479 


0.005U89 


(J. 8155177 




4.5000 


0.0994396 


1450. 1022034 


0.0000686 


0.8491784 




4.5500 


0.0204393 


803.9247055 


0.0000254 


0.8894575 




4.6000 


0.0490170 


1102.0391388 


0.0000445 


0.9939513 




4.6500 


0.0819387 


1413.5077515 


0.0000580 


0.9235846 




4.7000 


0.2198924 


2064.6186829 


0.0001065 


0.2474822 




71775100 


“O."0OT087(5 




0.0050071 


0. 0356658” 




4.8 01)0 


0.0292708 


968.6815491 


0.0000302 


0.8483709 




4.85UC 


0. 1724146 


2569.9429626 


0.0000671 


0.0044560 




‘».90CO 


C. 2382646 


2537.3323975 


0.0000939 


0.00154 76 




4.V500 


0. 1540685 


2235. 0189209 


0.0000689 


0.9586740 




•7.0 OOO 


0.1996510 


2663.2410270 


0.0000750 


0.9144567 





i(/t‘N IN FRACTION OF A CIRCLt 



92 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 

2. Library 

Naval Postgraduate School, Monterey, California 2 

3. Commander, Naval Air Systems Command 1 

Navy Department 

Washington, D. C. 20360 

4. Dr. E. S. Lamar 1 

Chief Scientist, Naval Air Systems Conmand 

Navy Department 
Washington, D. C. 20360 

5. Mr. G. L. Desmond 1 

Aerodynamics and Structures Technical Administrator 
Naval Air Systems Command 

Navy Department, Washington, D. C. 20360 

6. Commander D. M. Layton 1 

Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

7. Professor E. J. Andrews 1 

Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

8. Mr. Walter McNeil 1 

Flight Test Division, NASA Ames Research Center 
Moffett Field, California 94035 

9. Mr. William Crawford 1 

Computer Facility, NASA Ames Research Center 

Moffett Field, California 94035 

10. Lieutenant John K. Ready 1 

Carrier Suitability Branch, Flight Test Division 
Naval Air Test Center, Patuxent River, Maryland 20670 

11. Lieutenant Paul S. Norton 1 

Naval Test Pilot School 

Naval Air Test Center, Patuxent River, Maryland 20670 



93 



INITIAL DISTRIBUTION LIST (Continued) 

No. of Copies 



12. Profliesor A. E. Fuhs 1 

Department of Aeronautics 

Naval Postgraduate School, Monterey, California 93940 

13. Chairman, Department of Aeronautics 1 

Naval Postgraduate School, Monterey, California 93940 

14. Lieutenant William M. McDonald 1 

SMC #2855 

Naval Postgraduate School, Monterey, California 93940 



94 



UNCLASSIFIED 



Security Classification 



DOCUMENT CONTROL DATA - R&D 

(Smcucity cimMmlUcmtion of title, body of abatract and Indexing annotation muat be entered when the overall report le claeelfled) 



1 ORIGINATING ACTIVITY (Corporate author) 

Naval Poatgraduate School 
Monterey, California 93940 



2a. REPORT SECURI TY CLASSIFICATION 

UNCLASSIFIED 



2b. GROUP 



3. REPORT TITLE 



THE DETERMINATION OF THE DYNAMIC RESPONSE OF A SMALL SWEPT WING JET FIGHTER 
TO ATMOSPHERIC TURBULENCE USING THE POWER SPECTRUM METHOD OF ANALYSIS 



4. DESCRIPTIVE NOTES (T^pa oi report and Inclueive datea) 

Thesis 



5 AUTHORfS^ rLaar namo. fl ret name. Initial) 

NORTON, Paul S. 

LT, U.S. Navy 



6- REPO RT DATE 

Deceaber 1967 



7a. TOTAL NO. OF PAGES 



92 



7b. NO. OF REFS 

lA 



Sa. contract or grant no. 



b. PROJECT NO. 



9a. ORIGINATOR'S REPORT NUMBERfS^ 



96. OTHER REPORT NOfS^ (Any other numbera that may be aaaignad 
thla report) 



d. 



10. AVAILABILITY/LIMITATION NOTICES 




11. SUPPLEMENTARY NOTES 

None 



12. SPONSORING MILITARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 93940 



13 ABSTRACT 



The power spectral techniques of generalized harmonic analysis were 
applied to deteraiine the phase and amplitude of the transfer function between the 
vertical component of turbulence and the resultant pitch rate of a small swept 
wing jet fighter. The power spectrum noalysia was based on actual flight test 
data. In addition, the amplitude and {^se of the aircraft's transfer function 
are determined theoretically by two different methods and compared with the re- 
sults of the power spectrum analysis. It was established that the power spectrum 
Mthod, as applied to random turbulence, can be used to determine an aircraft's 
response to this random input on a continuous baa is and that power spectrum 
techniques have considerable potential value in the study of aircraft gust loading 



DD 



FORM 

1 JAN 64 



1473 



95 



UNCLASSIFIED 



Security Classification 



UNCLASSIFIED 



Security Classification 



\ 4 

KEY WORDS 


LINK A 


LINK B 


LINK C 


ROLE 


W T 


ROLE 


W T 


ROLE 


W T 


Response to Turbulence 
Power Spectrum Methods 
Harmonic Analysis of Turbulence 






• 






‘,4 



DD ,”r..1473 < > 96 UNCLASSIFIED 



; 1 i 0 



Security Classification 



A 3 t 4 0 9 



;*'r* *Tj 



